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ABSTRACT 

We perform a global fit to ~ 5 , 000 radial velocity and ^111, 000 magnitude measurements in 29 photometric 
bands covering 0.3 //m to 8.0/im distributed among 287 Galactic, LMC, and SMC Cepheids with P > 10 days. 
We assume that the Cepheid light curves and radial velocities are fully characterized by distance, reddening, 
and time-dependent radius and temperature variations. We construct phase curves of radius and temperature 
for periods between 10 and 100 days, which yield light curve templates for all our photometric bands and can 
be easily generalized to any additional band. With only 4 to 6 parameters per Cepheid, depending on the 
existence of velocity data and the amount of freedom in the distance, the models have typical rms light and 
velocity curve residuals of 0.05 mag and 3.5 km s"'. The model derives the mean Cepheid spectral energy 
distribution and its derivative with respect to temperature, which deviate from a black body in agreement with 
metal-line and molecular opacity effects. We determine a mean reddening law towards the Cepheids in our 
sample, which is not consistent with standard assumptions in either the optical or near-IR. Based on stellar 
atmosphere models we predict the biases in distance, reddening, and temperature determinations due to the 
metallicity and we quantify the metallicity signature expected for our fit residuals. The observed residuals as a 
function of wavelength show clear differences between the individual galaxies, which are compatible with these 
predictions. In particular, we find that metal-poor Cepheids appear hotter Finally, we provide a framework 
for optimally selecting filters that yield the smallest overall errors in Cepheid parameter determination, or filter 
combinations for suppressing or enhancing the metallicity effects on distance determinations. We make our 
templates publicly available. 

Subject headings: Distance scale — Stars: variables: Cepheids — Stars: abundances 



1. INTRODUCTION 

The classical Cepheids play a key role in the extragalac- 
tic distance determination, both as a direct way of ob- 
taining distances of nearby galaxies and as a calibrator 
for other methods ( e.g. the review by ^Freedman & Madore 
2pi0a, and Madore & Freedmanl 1199 1: Feast & Catc hpole 
1997; Freedman et al. 2001; Benedict et al. 2007; Rie ss et al l 
2009 a.b. ,2011a b). This is enabled by the relatively tight re- 
lation between the pulsational period of a Cepheid and its lu- 
minosity (the PL relation). Accurate distance estimates using 
Cepheids require increasingly better understandings of sys- 
tematic effects such as extinction, composition (metallicity) 
and blending. 

For example, the light from Cepheids is typically ex- 
tinguished by dust, which can have different properties in 
differ ent galaxies (e.g. jCardelli et al. 1989; Lanev & Stobi^ 
l993VFalco et al. l999; Misselt et al. 1999; Motta et al. 2002; 
Gordon etal. 2003; Indebetouw et al. 2005; Marshall et al. 
2006; Laney & Caldwell 2007). Cepheid observations must 
be corrected for extinction before any subsequent analysis. 
This usually makes use of an assumption that Cepheids with 
a given period have a specific color and that there is a 
known, universal extinction law. This holds for the Baade- 
Wesselink and related methods as well as purely photomet- 
ric analyses. To alleviate these shortcomings, recent ef- 
forts have shifted from the optical (Hubble Key Project) to 
the infrared wavelengths {H, I and Spitz er IRAQ, where 
the effects of reddening are reduced (e.g. iMcGonegal et al.l 



Tgs llGieren etani2q05t iPersson et al.ll2004t iFreedman et al 



2008^ iMadore et all l2009t iNgeow et all 120101; IStorm et al 
,201 laiiblK There are, however, differences in the near/mid- 
infra red extinction law for different sightlines in the G alaxy 
(e.g.l Roman-Zufiiga et al.ll2007HNishiyama et al.ll2009h . 

Cepheid properties must also change as a function of chem- 
ical composition, but the associated observational signatures 
have been a matter of considerable debate. While the gen- 
eral conclusion of observational studies is that metal-rich 
Cepheids are brighter than their met al-poor counterpart s (e.g. 
Stothers 1988; Kochanek 1997; Sak ai et al.ll2"0 04: Macr i et all 
i2006t; ,Valle et al.. ,2009^: ,Shappee & Stanekl ' 1201 L) . results 
from stellar pulsation models and some observation s yield 
the opposite depende n ce (e.g. Fiorentino et al. 2002, 2007t 

2008; Romaniello et all 
Estimates of the metal- 
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Marconi et al .l l2005t iBono et al. 
2008; Freedman & Madore 2011). 
licity correction vary by almos t an order o f magn itude, as 
can be seen in Figure 14 of iGerke etal] (1201 ll) . It is 
also now clear, at least in the extragalactic context, that 
uncertainties in the appropriate metallicities and metallic- 
ity gradients are nearly as important as the actual metalhc- 
ity dependence of Cepheid parameters (see iBresolinI 1201 It 
IGerke et al.1 120111; IShappee & Stanekl 120111 and references 
therein). Furthermore, studies find di f ferent trends i n dif- 
ferent filter s (Bono et al. 200l 120101; iNgeow et al.1 1201 ll 
IStorm etaTlbOl l b). Many of these problems arise because 
the effects of metallicity have filter-dependent degeneracies 
with reddening and distance, making it crucial to fully re- 
cover all covariances in the fina l results, as ernphasiz ed by 
Gould ( 1994), KochaniS Sl99% . iSasselov et al] Sl99% . and 
Riess et al.i (|201 lal) . As a result, the question of metallicity 



effects on the Cepheid distance scale remains a controversial 
issue. The problem of blending then adds further complica- 
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tions for more distant galaxies dMocheiska et alj|2000l) . 

In many circumstances, light curve templates for the indi- 
vidual observational bands are needed to accurately determine 
the mean magnitudes of extragalactic C epheids because of 
their sparse light curves. 'Stet sonI (11996!) constructed V- and 
/c-band templates for the LMC, SMC, and Galactic Cepheids 
with periods between 7 and 100 days by fitting Fourier series 
to the light curves. He assumed that the amplitudes and phases 
of individual modes vary continuously with p eriod, thus 
forming the well-known Hertzsprung sequence (Hertzsprung 
[1926). Hend ry et a l. (1999), Ngeow et al. (2003), Tanviret al. 
yOOS) and lYoachim et al.l (120091) expanded on these mod- 
els using Principal Component Analysis (PCA) techniques to 
build new templates for analysis of HST Cepheid data. All of 
these templates are limited in the sense that they are defined 
only for a small set of filters (typically V and Iq), and there is 
no well-defined means of shifting them to other, similar filters, 
let alon e to very different wavelengths. Freedman & Madorg 
(1201 Obi) show that all light curves are related by a linear trans- 
formation which provides a path towards addressing this prob- 
lem. 

In this work, we address these issues by constructing a 
global model of Cepheid light curves and radial velocity 
curves in the Galaxy, LMC, and SMC. We self-consistently 
determine distances, reddenings, radii, and temperatures of 
individual Cepheids along with a global reddening law, the 
mean SED of a Cepheid and its dependence on temperature, 
and the phase variations of the radii and temperatures as a 
function of period. We self-consistently determine the un- 
certainties in all parameters including all their covariances. 
Our model enables us to construct a light curve template for 
an arbitrary filter given just a single calculable parameter for 
the relative contribution of radius and temperature variations 
to that filter. We evaluate metallicity effects on Cepheid ob- 
servations first from a theoretical point of view, and then we 
search for metallicity effects in the estimated parameters and 
residuals of the fit. 

The approach to th e problem is closely related to th e Baade- 
Wesselink method ( Baadd 119261: IWesselinkI \194 &) and as- 
sumes that the magnitude of a Cepheid at a given time de- 
pends only its distance, reddening, and instantaneous radius 
and temperature. Rather than carrying out a study of individ- 
ual Cepheids, we simply fit a global model to all the available 
photometric and velocity data . A simpler variation on thi s 
idea was recently proposed by i Freedman & Madord (1201 Obi) . 
This paper is organized as follows. In Section |2] we outline 
our physical model, priors, dataset and the method of fitting. 
In Section [51 we describe results of the fit, discussing in turn 
the fit residuals, global quantities and individual properties of 
the Cepheids. In Section[4l we discuss our results and search 
for signs of additional physics, principally metallicity effects. 
We review our findings and outline future directions in Sec- 
tion [5] An Appendix outlines a general quantitative approach 
to selecting filters for Cepheid studies. 

2. MODEL AND DATA 

In this Section we describe our physical model of Cepheid 
light curves and radial velocity curves in detail (Sectio n [2. Il l 
along with its connection to traditional methods (Section lZSI i. 
The priors needed to address degeneracies in our master equa- 
tion, such as the zero point of extinction, are described in ^2.31 
The data that are used for the fit are described in ^2. 41 and the 
fitting method is outlined in Section[ 



2.1. Physical Model 

For each Cepheid, the radial velocity v(f) and magnitude 
OT,(f) in filter / at a time t are given as 

m,(f) = M,. + fi + ^iEiB - y) - 5 logf^) - 2.5 ftlogf^^V) 



- id;;(f) 

v(f) = v — , 

p at 



(2) 



where M, is the absolute magnitude of a "mean" Cepheid with 
radius R{) = 10 and temperature To = 5400 K, /i is the dis- 
tance modulus, £ii is the ratio of total to selective extinction in 
filter /, E{B-V) is the reddening, R{t) and T(t) are the radius 
and temperature of the Cepheid, v is the mean radial velocity, 
and p = 1.36 is a projection factor for co nverting observed 
radial velocities to p hotospheric ones re.g. lBurkiet"an[T982l: 
iNardetto et al.ll2004l) . Although the exact value of p and its 
dependence on the pulsational period is a matter of debate, 
these uncertainties largely cancel because we homogeneously 
anal yze the three galaxie s in our sample with the same value 
of p dStorm et al.ll201 lab . The coefficients /?,■ are the logarith- 
mic derivatives of the spectral energy distribution (SED) with 
respect to temperature at the reference temperature Tq, 



dlogFi 



d\ogT 



(3) 



where Fi is the flux in filter /. All logarithms in this paper are 
base 10. 

The radii R and temperatures T depend on the time t 
through the pulsational phase cf) = (t- to) /P, where to is the 
reference time and P is the pulsational period. We express the 
time-dependent parts of Equations ([T])-([2| as 



log (^^)=p-+A25p(0), 
log(^^)=r+A2^r(0), 



(4) 
(5) 



where p and t are the mean logarithmic radii and tempera- 
tures of a Cepheid with respect to the "mean" Cepheid with 
radius Ro and temperature To, is the dimensionless ampli- 
tude, and (5/3(0) and Sricj)) are the period and phase dependent 
changes in the radius and temperature. In the above defini- 
tions, quantities with subscript / are different for each filter, 
while /i, E{B - V), v, p, f, A^, P, and to are different for each 
Cepheid. Only the radius and temperature changes dp and 6t 
are functions of time t. 

Not all Cepheids in our sample have enough data to allow 
for independent estimates of p and f , but they can still con- 
tribute to the determination of global variables. In analogy to 
the normal period-luminosity (PL) relations, we assume that 
the mean radius and temperature of a Cepheid are functions 
of the pulsational period. 



{p(P))=ap + bplog{ — 
{fiP)) = af + br\og 



P 

lOd 
TOd 



(6a) 
(6b) 



where {p) and (f) are the average radii and temperatures at 
period P, and the coefficients Op, bp, af, and bf are estimated 
during the fit. These relations have corresponding widths, dp 
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and (Tr, which in principle could be estimated during the fit as 
well, but for simplicity we leave them fixed at ap = af = 0.02 
based on the differences between canoni cal and non-canonical 
theoretical pulsational models of Bon o et aT] (|T998, 2005). 
Equations (|6al)-(|6bl) also allow us to connect our approach 
to the traditional methods based on PL relations and to de- 
fine template light curves at any period or wavelength (Sec- 
tion |2]2]i. 

We model the time-dependent components as a truncated 
Fourier series of order A^^ , 

i5p = 51 ['^p-j cos(27r7(/)) + Spj &m{2-K}(f>)\ , (7) 

6t = '^ [Crj COSilnjcp) + SrJ sinilTTjcj))] , (8) 
where the coefficients are normalized in phase and amplitude 



such that c 



p.i 



1 and s 



; 0. The normalization of the tem- 



perature template with respect to the radial template is a part 
of the solution since the relative amplitudes of the radius and 
temperature changes are captured in the magnitude of Cr and 
Sr- We use a Cartesian representation of the Fourier coeffi- 
cients {Cpj,Spj) and icr.j,Sr.j) in order to avoid the coordi- 
nate singularities which occur in a "polar" representation (i.e. 
6p{(f)) oc a j cos[2tt j (j) + bj]) when the amplitude of a Fourier 
mode vanishes (a ,■ — ?• 0) and the phase is degenerate, as is seen 
in the iStetsonI (119961) Cepheid templates. In order to model 
the period dependence of the radius and temperature varia- 
tions, we specify Cp.y, Crj, Spj and Srj at Np periods chosen 
in such a way that the number of Cepheids in a given period 
bin is approximately constant (except for the longest periods). 
Values of the coefficients for a particular Cepheid's period are 
obtained by linear interpolation on the grid of templates in 
log P. We choose Nf = 20 and Np = 19. 

2.2. Connection to Traditional Methods 

In this Section we relate our approach to traditional PL stud- 
ies. Within our approach, the mean absolute magnitude M, in 
filter i of a Cepheid can be constructed from Equation ([T]l as 



M,=M,-5p-2.5Af, 



(9) 



where p and f are different for each Cepheid. This means that 
the average PL relation is defined by substituting the mean 
radii and temperatures from Equations (I6all-(l6b]i into Equa- 
tion (|9]l to obtain 

Mi = Mi-5{p{P))-2.5l3i{f{P)) = 

= M,- (5flp + 2 .5 A-flr ) - (5b p + 2.5 pibr) log ( 1 0) 

Thus our model leads to a PL relation with a zero point of 
M,-5ap-2.5/3,fl^ and a slope of -(5bp+2.5f3ibr)- For uncor- 
rected radius and temperature deviations (cp and af) from 
the mean trends (Eq. ||6l), the scatter in the PL is (5^(7^ + 

2.5^/3?cr^)'''^, so there is a strong correlatio n of wavelength 
(0i), s lope and scatter about the PL (see Mad ore & FreedmanI 
l201lh . Similarly, the mean template light curve for a filter / is 



mi(cj)) = -5Sp(<j))-2.5l3iST((t>), 



(11) 



which has a mean of zero ((m, ) = 0) and is scaled to an ampli- 
tude of = 1. iFreedman & Madore. (■2010b.) argue that light 



curves in one band can always be constructed as weighted 
sums of those in two other bands. This is true by construction 
for our models, where to produce band 3 from bands 1 and 2 
one chooses a scale factor x such that (3^ = x/Si +(1 -x)/32 so 
that m3(0) = jicmi((/)) + (l -x);«2((/'). In general, this can be an 
extrapolation since one can choose bands such that x > 1 or 
Ji:<0. 

2.3. Priors 

The parameters of the model are obtained by minimizing 
the master constraint H defined as 

H = x^ + S, (12) 

where is the sum of the squares of differences between the 
model and the observed magnitudes m°^'^ and radial velocities 



all data 



nObs 



E 



, ,obs 



(13) 



where a is the measurement error, and S includes contribu- 
tions from all priors, and the sums are over all stars, filters 
and measurements. 

Although we fit Equation dTST l to a huge dataset, the pa- 
rameters in Equation (fT3] l suffer from degeneracies without 
the addition of priors on the wavelength dependent vectors 
Mi, and /3, and the absolute distance and extinction scale. 
For example, we can simultaneously shift the extinction by 
EiB-Y) E{B-V) + l^E(B-y) and the absolute magnitude 
vector M, Mi- Mi^E{B-V^ while keeping m, the same. 
More generally, M, is well-determined up to adding compo- 
nents proportional to a constant (a change in the distance or 
radius scale), S^.i (a change in the extinction zero point), and 
A (a change in the temperature zero point). This is true also 
for S^i. Because we assume that /?,■ is the same for the mean 
and time-variable effects of temperature, it is not subject to 
the same degeneracies. 

In order to address these degeneracies, we first define priors 
for the vectors M,, /?,, and We add priors on M,, /3, and Sii 
of the form 



^M, - E 



(14) 



where Mf"" is obtained by convolving the 7}) = 5400 K, Rq = 
lORp,, log^= 1.5, [ M/H I = 0.0 "mean" Cepheid based on the 
ICastelli & Kuruc j (120041) model atmospheres with a tophat 
function defined by the central wavelength and width of each 
filter The sum is over all passbands. We choose cr^ =0.1 mag 
for all filters where we know the zero points. For the re- 
maining filter^ we estimate the conversion factor based on 
our data and then set ajj. = 0.2 mag. As we will see in Sec- 
tion 13.21 this assumption is unimportant. We constrain the 
vector f3i by evaluating Equation (O as a derivative in tem- 
perature about this reference model. The width of this prior 
is o-ft = 0.ipf°\ We fi x the coefficient /3v = 5.14 to its 
ICastelli & Kuruc3 (l2004l) prior value in order to prevent an 
overall shift of /?, due to a weak degeneracy with p. Simi- 
larly, we add a prior on where ^P"" is obtained from the 

' We were unable to find zero points for the CTIO JHK filters and the 
Washington filters CMTiTi. 
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ICardelli et all (Il989h galactic extinction law with =3.3 
and = + is held fixed so that E(B - V) has the 
standard interpretation. The remaining coefficients are con- 
strained by an "uncertainty" (t^_., = 0.05.^,!'"". The value 
My = 3.3 was chosen to agree with the Hub ble Key Project 
jMadore & Freedmanlll99UlFreedman et aljr2001.) . but is ac- 
tually an additional source of uncertainty that could be ex- 
plored. 

The second set of priors is on the mean radii and tempera- 
tures (Eqs. Il6al - ll6bl ). First, we add two priors 



all stars k 
all stars k 



Pk-im) 



fk-{f(P)) 



(15a) 
(15b) 



where (piP)) and (f(f )) are the mean radii and temperatures 
at a period P, which are given in Equations (l6ali-(l6b]i. and 
which are estimated as part of the fit. The widths of these rela- 
tions are fixed at (jp = <jf = 0.02 based on differences between 
the ca nonical and non-canonic al theoretical pulsational mod- 
els of Bono et all (Il998ll2005l) . These priors drive a Cepheid 
onto the mean period-radius and period-temperature relations 
if there are insufficient data to independently constrain its ra- 
dius and temperatures. The second set of priors is on the co- 
efficients of {p{P)) and (f (f )) in Equations. (I6all-(l6b]l. These 
have the form 



0.02 



bp- 



^pnor' 



0.01 



(16) 



along with a similar set of priors for and bf- The prior 
values of 



flr'' = 0.843, bT' = 0.655, 



prior 



-0.019, 



p 

.prior , 



-0.080, 



(17a) 
(17b) 



were set to match the "canonical" solar-metallicity, period- 
radius and effec tive temperature-luminosity relations of 
[Bono et al. (Il998h and [Bono et al. (2005). We derived the 
temperature priors assuming the standard relation of luminos- 
ity, radius and effective temperature L = AirR^aT'^. The prior 
widths of 0.02 and 0.01 were chosen to roughly correspond to 
the differences between the solar-metallicity "canonical" and 
"non-canonical" models of Bono et al. ( 1998, 2005), as an es- 
timate of the systematic uncertainties in the models. However, 
the specific choices for the prior values and their widths are 
of little consequence for the actual results. The zero points ap 
and Qf are well-constrained by the data. The slopes bp and bf 
are less constrained because of the limited number of long- 
period Cepheids and our restricted period range. For these 
coefficients, the width of the prior is somewhat important. 

The third set of priors set the distance and extinction scales. 
In order to fix the distance scale, we assume that the LMC 
Cepheids occupy a thin d isk with an inclinatio n of 30.7° and a 
position angle of 15 1 .0° (INikolaev et al.l2004 ) with a distance 
modulus at the center of 18.50 mag. The distances to the indi- 
vidual Cepheids in the SMC are free to vary, but we impose a 

prior that Mgl™ =18 .90 niag with a scatter of o-u,,, ^ = . 1 mag 
(e.g. iHilditch etaHIIOOl Keller & Woodlllool . The LMC 
distance scale cannot be fixed using a prior that assumes a dis- 
persion around some mean value like the one we use for the 
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Figure 1. Effect of the priors on the high-order Fourier coefficients and tem- 
plates at adjacent periods on the overall goodness of fit. The contours indicate 
the excess (denoted with numbers at each contour) over a fit without the 
priors 5ho and 5ap on the smoothness of the templates. The inset plot shows 
excess along the line of approximately equal contributions of 5ho and Sap 
to the X" (logo"ap = 1.81og(T[,o, dashed hne in the main plot). The black filled 
circle corresponds to the final choice of (Tho and (ti,o, which give a x^ excess 
compared to no smoothing of about 5%. 

SMC Cepheids because it leads to period-dependent residuals 
in distance moduli and reddenings which absorb any differ- 
ences between the requirements of the data and the priors on 
radii and temperatures. Essentially, the PL relation is dictated 
by the radius and temperature priors rather than the data. This 
problem would be solved by including a sample of Cepheids 
that truly lie at a common distance modulus. To fix the ex- 
tinction scale, we impose a prior that the LMC Cepheids have 
a me an extinction of {E{B-V)) = 0.147 mag (lUdalski et alJ 
Il999h with a width o'£(B-v) = 0.02 mag. 

The effects of most of these priors on the final results are 
weak, because most of the values are ultimately controlled 
by the data. The key exceptions are the mean extinction in 
the LMC, the distance to the center of the LMC, and setting 
= 3.3 and ^g-^v = 1- It is possible to recognize where 
the data dominate a prior by the final values and uncertain- 
ties - if the final results match the prior and the uncertainty is 
comparable to the prior width, then the data added no infor- 
mation, while if the values have shifted and the uncertainties 
are markedly smaller, then the data dominate. 

The last set of priors concern the templates themselves. In 
order to minimize unnecessary "oscillations" in the templates, 
we add a prior to keep the high-order terms of the Fourier 
series small 

^ho = -rE/(^L+4j+<;+4,,), (18) 

h° ./=3 

and to minimize differences between adjacent period bins we 
set 

^ A'f-1 



+ {Cr.j-Cr.j+lf + (Sr.j-Sr.j+lf] 



(19) 



The relative strengths of priors Sho and S^p and their relative 
weights with respect to other components of the H are deter- 
mined empirically following Figure [T] We vary the values of 
fJho and dap, and we record the change in the relative to a fit 
with negligible values for the weights. The final choice of Cho 
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and (Tho gives an equal weight to both priors and leads to an 
increase in the of the fit by about 5% compared to using no 
smoothing. We choose these particular values of (Xho and Cap 
to see some smoothing of the templates while only introduc- 
ing a small increase in compared to having no smoothing. 
The exact choice is somewhat subjective, but also has no im- 
portant consequences for the results. 

2.4. Data 

In order to obtain well-determined templates, the prop- 
erties of individual Cepheids and the global parameters of 
the model solution, we require a sample of Cepheids with a 
large quantity of photometric and radial velocity measure- 
ments. The data include the huge database of OGLE-III 
y and / measurements from Soszyriski et al. (2008; 2010), 
the m ajor da tabases of .Berdnikov (1999) and^orvnva et al. 
(119921 Il996l 119981). th e large sample of near-IR measure- 
ments by iPersson et al.' ('2004'), and t he Spitzer IRAC mea - 
sur ements o f Freedman et al. (2008]), iMarengo et al.l (120101) . 
and'Ngeo w & Kanbuil (12010). The references for the photo- 
metric and radial velocity data are given in Tables [T] and |2l 
respectively. 

These data include photometry in 29 filters covering the 
wavelength range from 0.3 /^m to 8.0 /im, as summarized 
in Table |3] In the optical, we include the standard John- 
son UBV , Cousins (^/)c and Johnson (Rl)i bands. We also 
use Hipparcos and Tycho photometry (bands Hp, Z?t, and 
Vt), photometry in the Washington system (C, M, T\, and 
Ti), and measurements in the Walraven photometric bands 
((VKt/LBy)w). The Washington system was designed to pro- 
vide metallicity and tempe rature estimates for G and K giants 
dWallerstein & Helferlll96 6). The Walraven system was de- 
signed for studies of early-type stars (Walraven & Walraven 
Il960h . We convert the Walraven data from their default logj^ 
scale to magnitudes by multiplying the data by a factor of 
2.5. In the near-IR we in clude the JHK filters of the SAAO 
(lGlasslll973 ) and CTIO ( Elias et al.l 1982) systems separately. 
Finally, we include the Spitzer IRAC [3.6], [4.5], [5.8], and 
[8.0] bands. 

Since we are mostly interested in extragalactic uses, we re- 
strict our dataset to fundamental-mode Cepheids with P > 
10 d. We do not in clude ultra-long period Cepheids with 
P > 100 days dBird e t al. 2009). As we proceeded with the 
fit, we purged the data set of obviously wrong measurements 
and database errors. We also removed Cepheids with too few 
data, binary systems, strongly blended systems, and Cepheids 
with obvious period changes. In the end, we have 177.314 
photometric and 5,031 radial velocity measurements for 287 
Cepheids in 29 photometric bands. 

2.5. Fitting Method 

Minimizing Equation ( fT2] l based on the physical model 
given by Equations ([T])-([8l) is somewhat similar to princi- 
pal component analysis (PCA). Specifically, we are trying 
to decompose the measurement ma in filter / for Cepheid 
k into an unknown set of vectors M,-, /?,■, 6p, and 6t 
while simultaneously obtaining the coefficients of the expan- 
sion fik, E(B-V)k, pk, fk, and A|. In a PCA, we would de- 
compose a data point m, as a sum m, = ajeij of coeffi- 
cients aj multiplied by orthonormal vectors ej, ^ij^ik = Sjk, 
simultaneously determining aj and ej. In this paper, we 
fit a measurement of a Cepheid k in filter / as a sum 
niik = Y^jOikjeij, where from Equation ([U it follows that e, = 



{Mi, 1,^, -5, -5(5p,-2.5/3, -2.5AM and a^. = {l,/x^,,£(B- 
V)k,pk,Al,fk,Al}. The basis vectors ej are not orthogonal 
and do not have unit norms in order to maintain their physi- 
cal meaning. The lack of orthogonality then leads to degen- 
eracies, which we must control by introducing the priors dis- 
cussed in Section |231 Mathematically, however, the similar- 
ity to PCA means that the model is well defined and can be 
solved by standard iterative or minimization methods to yield 
a unique solution given the data and priors. 

To this end, we have developed a versatile program that al- 
lows individual variables to be switched on or off, to reinitial- 
ize the physical variables to their prior values, and to alternate 
fitting all variables simultaneously with fitting just the prop- 
erties of individual Cepheids. We assign a weight to each 
measurement calculated as a maximum of the reported mea- 
surement error (if available) and a rms scatter of all measure- 
ments in the given filter from that particular data source. The 
fitting procedure itself can proceed in two ways. First, all 
variables are minimized independently in each iteration us- 
ing analytic first and second derivatives through the conju- 
gate gradient method (Press et al. 1992). This method is fast, 
but does not provide any explicit error estimates. The sec- 
ond option is to construct the full covariance matrix, which is 
then inverted using the Cholesky decomposition. This method 
is much slower, but provides error estimates that include all 
the correlations of the model as well as the measurement er- 
rors. Given any reasonable starting point, the fits are stable 
and well-behaved for all parameters except the phase refer- 
ence time to, where the surface is more complex and a 
manual intervention is sometimes necessary. While we have 
kept periods P and period derivatives P = fixed, they can be 
included without difficulty, but are not as stable because the 
surface in P and P is not smooth. 

3. RESULTS 

In the following Sections we present the detailed results of 
our model. In Sec tion [3. II we discuss residuals to the global 
fit. In Section lTSl we discuss the global parameters: the zero 
point Mi, te mper ature dependence /3, and extinction vector 
Section 13.31 examines the individual parameters of the 
Cepheids. 

3.1. Residuals of the Fit 

Figure |2] shows the data and the best-fit models for six 
Cepheids in the Galaxy, LMC and SMC. For the left col- 
umn, we have chosen stars with a large number of datapoints 
in many filters and a typical fit quality, while the right col- 
umn shows Cepheids with fits among the worst 5% for the 
given galaxy based on the photometric per degree of free- 
dom. We see that the model fits the data very well, and that 
even sparse data can strongly constrain the fits given the light 
curve structures required by the better sampled bands (e.g. 
HV 1543, in the lower left panel of Figure |2]l. The common 
problems are noisy data (HV 6320, in the lower right panel), 
inadequacy of the templ ate (SU Cru, upper right panel, see 
the discussion in Section l4~2] i and a small phase shift between 
data from different sources (HV 879, middle right panel). The 
remarkable point, however, is that this simple physical model 
of a radius plus a temperature template combined with a single 
amplitude produces good fits for objects that have photome- 
try in many filters across a broad wavelength range as well as 
their radial velocity measurements. 

Figure [3] shows the residuals of the fits as a function of 
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Figure 2. Examples of light curves and radial velocity curves for Cepheids in the thi'ee galaxies (Galaxy at top, LMC at center, SMC at bottom). The fits are 
shown with grey solid lines. The left column shows Cepheids with data in many photometric bands and median-quahty fits, while the right column shows the 
Cepheids with poor fits to their photometric data. 
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Galaxy 



<A-u> 



-0.100 ± 0.052. a = 3.54 




<m> = -c 



- <Afl> = 0.000 ± 0.001, a = 0.070 -- <&B> = 0.007 ± 0.003, ct = 0.108 -- <AS> = -0.009 ± 0.003, a = 0.094 - 



~» , J.,f i.'ir.'4.ia-.''^',re:»-3r^-^-^^ 



<AV> 



0.001 ± 0.000, a = 0.053 



<AR> 



-0.001 ± 0.001, CT = 0.044 



<A/> = 0.002 ± 0.001, a = 0.049 



LMC 

:Av> = -0.000 ± 0.173, CT = 3.1i 



<Al'> = ^0.001, ± 0.001, CT = 0.056 



<A/?> = 0.002 ± 0.002, CT = 0.051 



<A/> = 0.000 ± 0.000, CT = 0.029 



SMC 

<Av> = 0.000 ± 0.458, ct =r454 



<AV> = -0.000 ± 0.001, CT = 0.057 



<hR> = 0.012 ± 0.003, CT = 0.059 



<A/> = -0.001 ± 0.000, CT = 0.033 











- <AJ> = -0.008 ± 0.001, a = 0.041 - 


^ <AJ> = 0.005 ± 0.001, a = 0.047 ^ 


- <AJ> = 0.021 ± 0.004, a = 0.065 - 










; <AH> = 0.004 ± 0.001, CT = 0.037 ^ 


^ <AH> = -0.001 ± 0.001, CT = 0.047 ^ 


, <AH> = -0.009 ± 0.004, a = 0.070 










, <AK> = 0.000 ± 0.001, CT = 0.041 , 


- <AK> = -0.003 ± 0.001, CT = 0.049 - 


- <AK> = 0.008 ± 0.004, ct = 0.074 , 
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Figure 3. Re.siduals with respect to the model fits for radial velocities (top row; in km s^') and the most common filters (remaining panels). The red lines in each 
panel are averages in bins of A</< = 0.05 and the error bars are the uncertainties in these averages. Each panel gives the mean value of the residuals, its uncertainty, 
and the dispersion of the data around the mean. This dispersion is a combination of measurement errors and systematic residuals. 
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Figure 4. Residuals for the Spitzer IRAC bands. The format is the same as in Figure[3] but we do not show the phase-binned averages due to the small number 
of data points. For comparison to the dispersion of the residuals about the mean, we also report the average photometric uncertainty of the data. 

phase for the radial velocities and the most common opti- 
cal and near-IR bands. Figure |4] shows the residuals for the 
Spitzer IRAC bands. There are two separate ways to evaluate 
these residuals. First, for the mean properties of the Cepheids, 
we want the mean residuals to be consistent with zero given 
their uncertainties. This is generally true to high accuracy and 
it is reassuring that the average residuals are essentially zero 
for all filters with few systematic trends in the residuals as a 
function of phase. They are not exact ly ze ro, however, which 
will be an important point in Section |431 Second, we would 
like the dispersion of the residuals to be consistent with the 
measurement uncertainties. This is more difficult to evaluate, 
because not all the data sources include error estimates and we 
used as a weight the maximum of the actual measurement er- 
ror and the rms residual for each band from each data source. 
However, all measurements in the Spitzer IRAC bands have 
associated errors and therefore we can compare them to the 
fit residuals in Figure |4] We see that in the LMC the mean 
data error bars are similar to the scatter around the fit, sug- 
gesting that our model fits the data very well. In the SMC and 
especially in the Galaxy, there is a noticeable overall offset 
between the data and the model in some Spitzer bands, which 
causes formally a larger scatter around the mod el. T he origin 
of the offsets will be discussed further in Section l4!2] Further- 
more, we see in Figure[3]that the residuals are not entirely free 
of structure. The most obvious structure is that seen in the Ic- 
band LMC and SMC residuals at phases between (p = 0.4 and 
0.7. These phases correspond to the fast rise from minimum 
to maximum light. We defer the discussion of these patterns 
in the residuals to Section|4] 




Figure 5. The temperature and radius templates at the period grid points as 
a function of phase = -fo)/P. The templates are shifted with respect to 
each other for the sake of clarity. The period increases from the bottom to 
the top, and the actual period values are given to the right. The vertical line 
segments in both panels have unit length. 



3.2. 



Global Quantities 

Figure|5]shows the radius and temperature templates for the 
anchor points of our period grid. The values of the coefficients 
are given in Table |4] We make the templates publicly avail- 



able! including a code to calculate the average template light 
curve and radial velocity curve for any period within the pe- 
riod range and any filter along with the radius and temperature 
changes for that band. We clearly see the Hertzsprung pro- 
gression in the temperature templates by following the shifts 

'' ihttp : / /www, astronomy . Ohio- state ■ edu/ -pe jcha/ cepheids] 
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in the position of maximum and the bump on the rising branch 
to earlier phases at longer periods. The changes in the ra- 
dial templates are more subtle, and the main feature is that 
the minimum becomes more sharply peaked as the period in- 
creases. 

The magnitude zero-points M,- in Equation ([T]l represent 
the flux of the "mean" Cepheid with radius 7?o and temper- 
ature T(). Comparing the fitted values to a theoretical model 
of a Cepheid atmosphere allows us to judge the results. In 
Figure |6] we convert the values of M, to flu xes and com- 
pare them to a supergia nt model atmosphere (lKuruczlll979l: 
ICastelh & Kurucz(l2004l) with our reference To = 5400 K and 
^0 = lO^J© at lOpc assuming a fixed LMC distance modulus 
of 18.50 mag. We see t hat ou r fit gives fluxes consistent with 
the iCastefli & Kuruczl (120041) model. There appears to be a 
small systematic offset in the sense that our data are about 
0.08 mag too bright. Such a small shift can be caused by a 
number of reasons such as slightly incorrect flux calibrations, 
a small change of the LMC distance or a slightly different 
mean LMC extinction. The errors on M, presented in Table [3] 
are much smaller than the assumed width of the prior, indi- 
cating that Mi is well constrained by our data. We note that 
the temperature scale of our model is controlled by the prior 
on Mi so the mean period-temperature relation (Eq. |15b| i ad- 
justs to match this temperature scale. This calibration leads to 
Cepheid temperatures that are systematicall y higher than the 
period-temperature relation from the lBono et al.ldT998ll2005b 
pulsational models. If we removed the M,- prior and instead 
fixed the coefficients in Equation (I6ali-(l6b]i. we would get per- 
fect agreement of the mean temperatures with their prior val- 
ues (Eq. Ill7bl ) at the cost of poor match to the spectral energy 
distributions (M,) of the lCastelli & Kuru cz (2004) model for 
that temperature. Clearly, having a theoretical stellar atmo- 
sphere model that is known to be correct would allow us to 
use a stronger prior to constrain our model and remove most 
of these ambiguities. 

As discussed in Section [23] (Eq. [S]), the temperature co- 
efficients /3i are the logarithmic derivatives of the flux in any 
band with respect to temperature at temperature Tq. Figure |7] 
compares the fitted values of to those for a 5400 K black 
body and to the prior based on the ICastelli & Kuruc d (l2004t) 
atmosphere models with \ogg= 1 .5 and at To = 5400 K. Recall 
that = 5 . 1 4 is fixed to its prior value based on the model at- 
mospheres. Our results clearly deviate from the black body in 
the sense predicted by the model atmospheres. At short wave- 
lengths, the Pi are larger than the black-body model because 
of the effects of metal opacities, and the model may even re- 
solve some of the expected spectral features. In the infrared, 
we clearly see the CO band head at ^ 5 fim. We obtained es- 
sentially the same result with a black body prior on /?,, which 
means that the weakly imposed prior on /3, is simply over- 
whelmed by the statistical power of the data - the values of /?,■ 
are robustly determined by the color changes of the Cepheids 
during the pulsational cycle. There are clear systematic dif- 
ferences from the Castelli & Kur ucz (2004) model, especially 
in the infrared, which are l ikely due to real p r obleni s in the 
molecular o pacities of the ICastelh & Kuruc d (120041) model 
dPremaux e t al. 2006). 

In Figure |8] we compare the final extinction curve to the 
ICardelli et al. (1989) .Siv = 3.3 extinction curve used as a 
prior. We focus on the shape of the extinction curve Mil^v 
because we fixed Siy = 3.3 in the models. In Figure [8] one 
should also focus on the narrower band pass filters because 



we have not convolved the extinction law with the average 
Cepheid spectrum and filter bandpass. For example, the very 
broad Hp filter has anomalously low presumably be- 
cause it is effectively a much redder band pass in the pres- 
ence of significant extinction. An advantage of the model 
S^i is that they "correctly" include all these band pass av- 
erages. Figure |8] shows that between 0.5 and 1.5/im"', the 
results fall below the prior, while for A"' > 2/im"' they are 
above. This suggests that the "mean" extinction law of our 
ph otometric datas et is different from the canonical = 
3.3 ICardel het al.1 (IT989.) law, fafling between the = 3.3 
ICardelli et all (Il989hlaw and the empirical LMC and SMC 
curves of Gordon et al.l (I2003h . The S^i values are controlled 
by the data except for the Spitzer bands, where t he values and 
uncertainties closely follow the 'Cardelli et al. (1989) prior. 
That the best fit extinction curve is not simply a Cardelli et al] 
(1989) = 3.3 law should not be a sur prise. There is n o 
unique extinction law even in the Galaxy (iBerrv et al.ll201 Ih . 
and the L MC and SMC extinction laws show furt her differ- 
ences (e.g.l Misselt et al.lll999t iGordon et al.ll2003h . Here we 
have determined the best average extinction law for Cepheids 
in these three galaxies. We can att empt to param eterize our 
results within the framework of the ICardelli et al.l (.198 9) ex- 
tinction models. Becau s e we keep fixed, we perform a 
fit of the ICardelli et al.l (119891) extinction law evaluated for 
Acenter of filter i to S^ij Siy with ^y as a free parameter We 
find that the best fit value is ^y = 3.127 ±0.002, but most 
of the Mil Sly values still show significant offsets of several 
percent from the best fit value. Interestingly, the largest dif- 
ferences are for the near-IR bands, where /Tsaao and Actio are 
shifted by about 47% and 31% from their respective best-fit 
'Cardelh et alj (Il989h values, as can be seen also in the inset 
of Figure |8] We also fit the Cepheids data using M.y =3.1. 
The change in is negligible and the ,^,/^y curve tilts to 
a slightly steeper slope to accommodate = ^y + 1, but 
it is still discernibly different from the SiijSiy ratios of the 
^y = 3 . 1 ICardeUi et al.l (fT989l) reddening law. 

After establishing that /?,■ and represent the physical 
quantities intended in Equation ([U, we can assess whether 
temperature and extinction can be determined independently. 
In Figure |9] we show the ratio of /3; to as a function of 
inverse wavelength. Because changes in temperature and red- 
dening can have similar effects on the observed colors for a 
limited wavelength range, a change in f can be mimicked by 
a change in E{B-Y^ if I5i/^%i is constant over that range. We 
see from Figure |9] that the ratio Pi/^%i has an approximately 
parabolic shape in A"' , and therefore we should be able to sep- 
arate the effects of reddening and temperature quite robustly if 
there is enough wavelength range. Determining both temper- 
ature and extinction requires a minimum of three filters (two 
colors) and the efficacy of any choice can be assessed by con- 
necting the longest and shortest wavelength bands by a line in 
Figure|9l and then examining the distance of the middle band 
from the line. Temperature and extinction degeneracies are 
minimized by maximizing that distance. We generalize this 
procedure in Appendix IaT In assessing Figure^ it is impor- 
tant to recognize the large vertical scale. The curvature be- 
tween filter wavelengths need only be large compared to the 
photometric errors (< 0. 1 mag). For the same errors, however, 
a larger wavelength baseline is almost always better 

We can check for possible temperature/extinction correla- 
tions by looking at the individual temperatures and extinctions 
shown in Figure[TO] In the LMC, there is a group of stars with 
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Figure 6. The luminosity of the "mean" Cepheid with /?o = lO/?0 and To = 5400 K located at lOpc constructed from the values of M; and converted to \F\ using 
the LMC calibration distance of 18.50 mag. The horizontal error bars show the width of each filter and the filter names are given at the top. The uncertainties 
on Mi are shown, but they are generally smaller than the size of the symbol. For comparison, we also show an error bar corresponding to the width of the 
prior on M,. The grey lin e shows the flux (\F\ at lOpc) of a supergiant model atmosphere with To = 5400 K, Rq = WRq, and \ogg = 1.5 interpolated from the 
ICastelli & Kuruca 420041) grid of atmosphere models. 



only V and Ic photometry, which does not allow for indepen- 
dent extinction and temperature determination. These stars 
concentrate along the prior {E(B-V)) = 0.147 mag. There 
is also a tail of LMC and SMC stars with both high f and 
E{B-V). These stars lack photometry in bands bluer than 
V, so part of the problem may lie in degeneracies created by 
small systematic effects in the photometry calibration and the 
extinction law. For the most extreme cases the model pre- 
dicts U -band magnitudes at maximum brighter than in the B 
band, which cannot be true. We consider the temperatures 
and extinctions well separated for most of the stars, because 
Cepheids in the Galaxy have much higher extinctions than 
those in the LMC and SMC but noticeably lower tempera- 
tures. Metallicity differences also introduce apparent shifts in 
temperature and extin ction. Arrows in Figure[TOl based on the 
results of Section l4n show the effect of increasing the metal- 
licity by 0.5 dex, which is approximately the spread in LMC 
Cepheid metallicities and the typical dif ference between LMC 
and SMC Cepheids dRomaniello et alj |2008). While the am- 
plitude and sense of the shifts are consistent with metallicity 
effects, none of the outliers has a measured metallicity to ver- 
ify this origin. 

Some choices for our priors are unimportant as the data so 
strongly constrain the problem that the same final results will 
be obtained even with very different choices for the priors. 
For example, using a black body prior for /3, instead of the 
estimate based on the ICastelh & Kurucd ( 120041) model atmo- 



spheres changes the results by a negligible amount. The de- 
partures from a black body due to metal line blanketing and 
molecular opacities are a robust result, as can be seen by com- 
paring the fit errors on the global parameters with the prior 
widths in Table [3] (prior widths are 5% for ^„ 10% for /?,, 
and 0.1 or 0.2 mag for M,, while the typical fit uncertainties 
are typically 1% for Mi, 0.7% for /3, and 0.01 to 0.02 mag for 
Mi, respectively). Our choice of priors is also overdetermined. 
For example, the prior on M, essentially drives o ur tempera- 
ture scale to be the same as in ICastelli & Kuruca (2004). but 
the prior o n f (Eg. [161) drives it (weakly) to the temperature 
scale of thel Bono et alj ( l2005h pulsational models. Similarly, 
radial velocities and multi-filter photometry implicitly define 
the distance to the LMC through the Baade-Wesselink method 
and so in theory no prior on the LMC distance would be nec- 
essary. We discuss these issues further in the Conclusions 
(Section|5l). 

Finally, we check our assumption that the LMC Cepheids 
are distributed in a thin disk with the inclination and posi- 
tion angle given by Nikolaev et alJ (2004) at an LMC center 
distance of /i = 18. 50 mag. In Figure [TT] shows the /c-band 
residuals for the LMC Cepheids as a function of their posi- 
tion angle relative to the center of the LMC. An incorrect in- 
clination or position angle of the disk would produce excess 
residuals at some position angles, but we see very little varia- 
tion in the scatter of the Cepheid mean residuals as a function 
of the position angle. We note, however, that modest inconsis- 
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Figure?. The temperature coefficients ft (Eq. (3]) as a function of wavelength. The horizontal error bai's show the widths of the filter passbands and the 
uncertainties in the ft are smaller than the size of symbols. The solid grey line is the prior on /3, based on the Castelli & Kurucz (2004) atmosphere models 
with To = 5400 K and logg = 1.5. The widths of the prior are illustrated by the error bars. The dashed grey fine is the estimate of ft for a black body with 
To = 5400 K. The biggest deviations from the black-body agree well with the Castelli & Kurucz (2004) model. There are, however, systematic disagreements 
with the Castelli & Kurucz 12004 ) model in the infrared that are probably caused by problems in the molecular opacities. 



tencies in distance can be partially absorbed by changes in the 
mean radius, especially for Cepheids without radial velocity 
measurements. 

3.3. Individual Properties of Cepheids 

In this Section we discuss our results on the individual prop- 
erties of Cepheids, which are given in Tables |5] |6] and |2l 
First, we compare our distances for the Galactic Cepheids 
with measured parallaxes to check our assumption that the 
LMC distance modulus is 18.50 mag. The overlap between 
our sample and the sample of Benedict et al. (2007), who 
measured Cepheid parallaxes using the Fine Guidance Sen- 
sor on the Hubble Space Telescope, is only 2 stars, £ Car and 
C, Gem, because o f our restriction toP ^ 10 days. The parallax 
measurements of [Benedict et al.l (l2007h give /i = 8.48 ± 0.22 
and 7.78 ± 0.14 for these two stars, respectively. We obtain 

= 8.29 ± 0.01 and 7.89 ± 0.03 for ^ Car and C Gem, respec- 
tively, given our assumed LMC distance of = 1 8 .50 mag. We 
note that the uncertainties in our distances include all the un- 
certainties in the model (ie. including /?,, etc.) and are not 
simply the uncertainties in the distance for fixed model param- 
eters. For these two stars, the distances are in good agreement 
w ith their p arallaxes. 

IStorm et al. (201 la b) determine Baade-Wesselink/I R dis- 
tances and extinctions (some from lPersson et al]|2004j) for a 
sample of Galactic and LMC Cepheids. For the 21 LMC 
stars we have in common, our distance moduli are higher 
by 0.054 ± 0.043 mag, where the error is the uncertainty in 



the mean, consistent with their lower mean LMC distance 
modulus of ^uAC = 18.45 ±0.04 mag compared to our fixed 
mean value of /ilmc = 18. 50 mag. Our LMC extinction esti- 
mates are also systematically higher by AE{B-V) = 0.036 ± 
0.020 mag. In the Galaxy, we have 30 stars in common, 
and in this case our distance moduli are shorter by an aver- 
age of 0.140 ±0.058 mag. Given that our LMC distance is 
about 0.05 mag longer, our Galacti c dista n ce scal e is about 
0.19 mag shorter than that of Stor m et al.1 (1201 lal) . For the 
Galactic Cepheids, we find slightly lower average extinction 
with AE(B-V) = -0.015 ±0.009mag. The LMC extinction 
offset could be reconciled by setting the mean LMC extinc- 
tion to (E(B-V) ) - 0.10 rather than the {E(B-V)) = 0.147 
we adopted from lUdalski et al.l (119991) . However, lowering the 
mean LMC {E{B-V)) to 0.10 mag would create a noticeable 
overall shift of the PL relation. 

Fi gure [T2] shows the distribution of the Cepheids in radius 
(Eq. Ill5al ). The overall distribution of the Cepheids tracks the 
prior well, but most Cepheids have errors in their radii that are 
smaller than the width of the prior, which indicates that the 
radius is constrained by the data (or other priors) more than 
by the prior on the radius. The larger scatter in the Galactic 
Cepheids is caused by a much greater spread in distances and 
reddenings compared to the LMC and SMC. A potential ori- 
gin for the large positive outliers is blending, where additional 
light from an unresolved companion causes the amplitude to 
be smaller and the Cepheid to be brighter. During the pro- 
cess of cleaning our dataset we removed one obvious outlier 
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Figure 8. The ratio ,^, /.^v/ of the total to the selective extinction Mi in filter i relative to .^v as a function of inverse wavelength. The horizontal error bars show 
the widths of the filter passbands, and the errors in the estimates of are smaller than the size of the symbol. The My =3.3 Cardelli et al. ( 1989) prior is shown 
by the solid grey line. The widths of the prior at the filter wavelengths are indicated by the error bars at the bottom of the plot. The dashed gray line shows an 
Mv = 2.5 Caidelli et al. ( 1989) reddening curve, and the gray dotted (dash-dotted) lines show the LMC (SMC) reddening curve of Gordon et al. ( 2003). The inset 
shows the infrared region in more detail. The meaning of the grey solid and dashed lines is the same as in the bigger plot. The red dotted and blue dash-dotted 
lines show the infrared extinction curves of Nishiyama et al. (2009) and Roman-Zuiiiga et al. (2007), respectively. The curves in the inset were normalized to our 
fitted value for .^^/Csaao ■ ^'^^ IRAC values of .i^, are dominated by prior, while the JHK values are not. 
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Figure 9. The ratio fii/l^i of the temperature coefficient /3, to the reddening 
factor .^i . Horizontal error bars show the widths of the filter passbands and 
the uncertainties in the ratios are generally smaller than the size of the sym- 
bol. The gray solid lines are the predictions for the Castelli & Kurucz (2004) 
and black-body expressions for /3, combined with an ,'%v =3.3 Cardelli et al. 
(1989) law for .5?, . The parabolic shape indicates that our span of wavelengths 
is large enough to ensure that the temperature and extinction estimates are not 
mutually degenerate. 



(HV 2326) created by blending (ISoszvnski et al.ll2008h . Our 
fit also provides the updated values of coefficients Op and bp, 
which describe how the mean radii change with the period. 
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Figure 10. Reddenings E(B — V) as a function of temperature f relative to 
the mean period-temperature relation (t) for the Galaxy (red open squares), 
LMC (green filled circles), and SMC (blue open triangles). The horizontal 
grey dashed hne is the LMC reddening prior {E(B — V)) = 0.147 mag with 
a width of 0.02 mag and the vertical dotted lines show the width of the (r) 
prior The grey arrows show the expected change of f and E{B — V) as a 
result of a metallicity increase of 0.5 dex due to changes in the atmospheric 
metallicity at fixed temperature (lower arrow) and changes in the physical 
properties of the star a t fixed period based on pulsational models (upper ar- 
row). See Section l4~n for details. 

We find 



(p) = (0.777 ± 0.002) + (0.698 ± 0.006) ^og[—j. (20a) 
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Figure 11. The /c-band residuals of the LMC Cepheids (dots) as a function 
of the position angle of the Cepheid relative to the center of the LMC. The red 
circles show the mean residuals for each Cepheid where the error bar is the 
uncertainty in the mean. The solid black line shows the expected signal in the 
residuals for Cepheids at the observed mean distance from the LMC center 
of 1 ° if we did not correct for the LMC tilt. Vertical solid and dashed lines 
mark the position angles of the line of nodes and the perpendicular direction, 
respectively. 
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Figure 12. Radii of individual Cepheids as a function of period for the 
Galaxy (red open squares), LMC (green filled circles), and SMC (blue open 
triangles). The grey solid line is the final mean period-radius relation {p(P)) 
of our dataset (Eq. I20al ) with the width indicated by the grey dotted lines. 
The g rey dashed line is the pri or on (p(P)) based on the .Bono et al.. 1 19981 
120051) pulsational models (Eq. I17al ). 

Compared to the prior, we find that the zero point is lower by 
about 0.07 dex, about three times the width of the prior on Op, 
and the output error in the zero point is much smaller than the 
width of the prior (0.02), which suggests that the data con- 
strain it well. On the other hand, the value and uncertainty in 
the slope are more similar to the prior, which suggests that the 
prior is important for the slope estimates. Our limited period 
range and the small number of long-period Cepheids are not 
particularly well-suited for the determination of slopes. 

Figure [O] shows the distribution in mean temperature. The 
final zero point and slope of the period-temperature relation 
are very different from the theoretical prior (Eq. Ill7bl ) based 
on the Bono et al. ( 1998, 2005) pulsational models. The zero 
point offset generally stems from the fact that our tempera- 
ture scale is constrained by the prior on M,, which does not 
exactly match the theoretical pulsational model of (f(P)), as 
was discussed in Section [321 We again see that the tempera- 
ture error bars are frequently far smaller than the width of the 
prior and so must be tightly constrained by the data. Figure [T3] 
also shows that the SMC Cepheids seem to be systematically 
hotter than their LMC and Galactic counterparts. We will con- 
sider the differences between galaxies in more detail in Sec- 
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Figure 13. Same as Figure |T2] but for temperature. The grey solid line show s 
the final mean period-temperature relation {f{P)) of our dataset (Eq. l20bl ) 
with its width denoted by the gray dotted l ines. The gr ey dashed line is the 
prior on ( f(P)) based on the IBono et alj <1998l |20()5I) pulsational models 
(Eq. [173). 

tion[4l Our fit also provides the updated values of coefficients 
Of and bf, which describe how the mean temperatures change 
with the period. We find a final period-temperature relation 
of 



(f) =(0.031 ±0.002) -(0.061 ±0.005) log l^—j . (20b) 

These values for and bf diffe r noticeably from the prior 
based on the iBono et all ( 119981 [2005) pulsational models 
given in Equation ( fTTI l in the sense that we find a higher zero 
point and a shallower slope. The uncertainties in both quan- 
tities are smaller than the width of the prior, which suggests 
that the data c onstr a in the se parameters well. 

Equations ( I20al )-( l20bb provide priors for the mean radii 
and temperatures for stars which do not have enough data 
to constrain these parameters independently and also define 
the mean PL and template light curves as outlined in Sec- 
tion l2.2l By combining the radii a nd te mperatures of individ- 
ual Cepheids as shown in Section IZ2I we can construct a PL 
relation as defined in Equation JTOt . We show the V-band PL 
relation in Figu re [T4l along with the the oretical PL relation 
impli ed by the Bono et all (119981 l2005h pulsational models 
(eqs. lll5aH15bl ). whic h has slope of -2.25. W e also show the 
y-band PL relation of lUdalski et al.l (119991) with a slope of 
-2.76 that was obtained by dereddening the LMC Cepheids 
using nearb y red clump stars. We se e that despite having the 
priors from iBono et al.l (119981 12005 ) theoretical models, the 
data control the fit in the sense that the final PL relation of 
the stars in our sample is now more similar to the empirical 
PL oflJdalski et al. ( 1999). The theoretical period-radius and 
period-temperature relations from Equation (ITtI i predict a PL 
relation t hat is incompatible w ith both our data and the earlier 
results of lUdalski et al.l ( 119991) . It is important to understand 
that the actual mean magnitudes of each Cepheid in any band 
are essentially equal to the mean of the observations, as seen 
in Figure [3] and that any differences in Figure [l4l arise only 
from different assumptions about the mean extinctions and 
distan ces. To pu t this on more quantitative grounds. Equa- 
tions (I20ab - (l20bb together with Equation ( ITOl i imply a PL re- 
lation 



Mv = (-4. 1 30 ± 0.024) - (2.703 ± 0.070) log 



lOd 



(21) 



where My is the absolute magnitude in the V band. The er- 
rors of the slope and zero point consistently include the un- 
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Figure 14. The V-band period-luminosity relation for t he C epheids in our sample. The meaning of the symbols is the same as in Figure mi The dot-dashed gray 
line is the LMC PL relation of Udalski et al . ( 1999) (Eq. 1221 ) assuming /^lmc = 18.50 mag and corrected for extinction based on red clum p star s . The dashed line 
is the PL relation defined by Equation UOt constructed from the priors on the mean period-radius and period-temperature relatio ns (Eg . Il7al - I17bl ). The solid 
grey line show s the "mean" P L relation constructed from the final period-radius and period-temperature relations in Equations )20a) - t20bl . Our output PL is 
consistent with 'Udalski et al.'fi999) even though the priors are not, which suggests that the Bono et al. 1 1998, 2005) priors are inconsistent, but are overwhelmed 
by the statistical power of the data. The arrows show the expected change of My as a result of a metallicity increase of 1 .0 dex due to changes in the atmospheric 
metallicity at fixed temperature (left aiTow) and changes in the physical properties of the star at fixed period based on pulsational models (right an'ow). 



certainties in My, ap, bp, Of, and bf and all their mutual 
covariances. The standard deviation of the Cepheids about 
this relation is 0.46 mag. Longer wavelengths show less scat- 
ter, with a standard deviation of 0.33 mag at Iq band and 
only 0.22 m ag for the IRAC [3.6] band. For comparison, 
lUdalski eta l. ( 1999) give a standard deviation for the V-band 
PL relations of 0.16 and 0.26 mag for the LMC and SMC, 
respect ively, and 0.11 and 0.2 1 mag for the Ic band, respec- 
tively. iScowcroft et al.l (1201 Ih gives a scatter of 0.14 mag in 
the Spitzer [3.6] band for the P > 10 days LMC Cepheids 
based on phase-resolved Spitzer photometry. Thus, the scat- 
ters about our PL relations are larger, but we are far more 
tightly constrained due to the large number of bands and we 
expect some additional contributions to the scatter from miss- 
ing physics because our sample also includes stars from the 
LMC, SMC, and the Galaxy, as we will discuss in Section H] 
We point out that given the limited period range of our sam- 
ple and the lack of long-period Cepheids, the slopes of rela- 
tions are poorly constrained. Adding short-period Cepheids or 
long-period stars in external galaxies would better constrain 
slopes of our period-radius, period-temperature and conse- 
quently period-luminosity relations. 

The absolute V-band magnitude of a Cepheid with P = 
lOdays (Eq. 1211 ) is approximately -4.13 ±0.02 mag, which 
is somewhat fainter than the -4.22 ± 0.04 found from the 



lUdalski et al] (Il999h V-band PL relation of 

Mv. = (-1.458 ±0.021)-(2.760±0.031)log^^^ . (22) 

As we experimented with our models, we would find small 
offsets between the zero point of our PL relation and the 
OGLE PL. These offsets are very sensitive to the definition 
of the E(B - V) extinction zero point in the LMC. For exam- 
ple, shifting the LMC extinction zero point by AE{B - V) = 
0.03 mag from our default {E(B-V)) = 0.147 mag offsets the 
PL zero point by ^vAE(B-V) = 0.10 mag. We also see in 
Figure [14] that there are potential differences between the in- 
dividual galaxies such that a PL fit only to the LMC Cepheids 
will be in better agreement with Equation ( |22] |. We discuss 
this in more detail in Section l43] 

FigurefTSlshows the amplitude (eqs. ||3S) as a function 
of period. There is a trend with the period: amplitude 
increases with period up to P ^ 30 d and then either declines 
or becomes constant at higher periods. Interestingly, for a 
fixed P, the distribution of appears to be skewed towards 
lower values of A^. In the case of LMC and SMC Cepheids, 
blending can decrease the amplitude somewhat, but the 
same behavior is seen is seen also in the Galactic Cepheids. 
We also do not see any correlation between A^ and p- (p). 
This suggests that a fraction of Cepheids have intrinsically 
smaller amplitude, presumably due to their position in the 



Global model for Cepheids 



15 



□ Galaxy 
• LUC 

A SMC 



30 40 
Period [days] 



50 60 70 ao 90100 



Figure 15. Amplitude as a function of period. Tlie meaning of tlie sym- 
bols is the same as in Figure 1121 

instability strip (Szabo et all 120071) . There are small differ- 
ences between the individual galaxies, which we discuss in 
Section 14.41 Our findings are in agreement with previous 
result s (e.g. [S andage & Tammann 1971; van Genderen 
1978t iBer dnikov & Ivanoy. ,1986: Paczvnski & Pindqr, 
2000t [Saiid age et a l.l I2004t iKlagyivik & Szabados 200a 

2om. 



iBer d 
iSand a 
Szabados & Klagyivi] 



4. DISCUSSION 

Our physical model, which includes only radius and tem- 
perature, does not include all relevant Cepheid physics - in 
particular it does not include composition. In this Section 
we discuss the evidence for additional physics that was not 
included in Equations ([T]i, (|2|i and We are interested in 
differences between the Cepheids i n the three galax i es be- 
cause of their different metallicities. iRomaniello et alj (120081) 
find mean spectral abundances of [Fe/H] ~ 0.00, -0.33 and 
-0.75 dex for the Galactic, LMC and SMC Cepheids with 
typical ranges of -0.1 8 < [Fe/H] < 0.25, -0.62 < [Fe/H] < 
-0. 10, and -0.87 < [Fe/H] < -0.63, respectively. It is impor- 
tant to remember that any effect of composition that can be 
mimicked by a parameter in the model has been! For example, 
a zero point difference M,- between the galaxies due to metal- 
licity effects is automatically compensated for, because we fit 
for individual distances /i for the SMC and Galactic Cepheids. 
Similarly, differences in color are absorbed into reddening and 
temperature before leaving any trace in the residuals. We start 
by investigating these issues using theoretical model atmo- 
spheres in Section l4n Bearing these issues in mind, we can 
search for extra p hysi cs on three levels. First, in the model 
residuals (Section l472] l. second, in global statistical properties 
like distances and PL relations (Section 14.3b . and, third, in 
the distributions of the properties of t he individual Cepheids 
(Section [4.4b . Finally, in Section |43] we evaluate the signif- 
icance of second-order terms and the first-order metallicity 
correction relative to the first-order temperature term 2.5/3f 
in Equation ([T]). 

4.1. Theoretical Predictions for the Effects of Metal Content 

Composition affects both the stellar atmospheres (chang- 
ing Eq. lHI) and the mean properties at fixed period (changing 
Eq. 0). The observed properties will be a combination of 
both effects. We can get a sense of metallicity effects on the 
mean Cepheid properties by examining model period-radius 
and period-ternperatu re relations at different metallicities. In 
the iBono et aTl (11998 ) models the radius at fixed period in- 
creases as the metallicity decreases, by 6% for Cepheids with 
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Figure 16. Magn i tude differences as a function of wavelength between 
ICastelli & Kuruczl <2004h stellar atmospheres with [M/H] = -0.5 and 0.0 
at fixed temperature To = 5400 K, radius Rq = IORq, and surface gra vity 
logg =1.5 (black line). The red line is the best fit our model of Eq. i2i\ 
can make to the effects of metallicity. 

P = 10 days when the metal content is decreased from Galactic 
to SMC me t allicit y. Combining the period-radius relation of 
iBono et aTl (Il998h w ith the temperature-luminosity relation 
of iBono et al.l (120051) . we find that the temperature increases 
by 5% for Cepheids with P = 10 days when the metal content 
is decreased from the Galactic metallicity to that of the SMC. 
In combination, t he luminosity increase s by ~ 35% at the 
lower metallicity. iFiorentino et al.l (120021) and iMarconi et al.l 
( 120051) also investigated the effects of metallicity on the mean 
properties using pulsational models, finding that increasing 
the metal content moves the instability strip to lower temper- 
atures. These effects will be absorbed into shifts in f and fi/p. 

Metallicity also modifies the SED of the stellar atmosphere 
at fixed T, L and R by changing the strengths of spectral lines 
and the a mount of line blanketing. Figure [16] compares the 
flux of a ICastelli & Kuruc j (120041) model atmosphere with 
[M/H] = -0.5 to one with solar metallicity at fixed temper- 
ature and radius (To = 5400 K and R = 10 Rq). As expected, 
the lower metallicity star is brighter blueward of ^ 0.45 /im, 
due to the reduced effects of metal-line blanketing, and in the 
molecular band heads at ^ 2.5 /im and ^ 5 /im. Now let us 
assume that we are fitting photometric measurements of stars 
with different metallicities using the model in Equation ([T]i, 
which does not explicitly take metallicity into account. This 
leads to shifts in distance, reddening, and temperature, be- 
cause these quantities will try to absorb as much of the vari- 
ance in the appearance of the stellar atmosphere as possible. 
To quantify this effect, we fit the magnitude difference Am 
between model atmospheres with different metallicities (Fig- 
urefTSIl with a model based on Equation ([T]), 



Am= Afi+^AE(B-V)-2.5(3Af. 



(23) 



For the purposes of this discussion, changes in /i are degener- 
ate with changes in p, and we thus include o nly A/i. We also 
assume the = 3.3 ICardelh et all ([T989h reddening curve 
for and the (3 computed from the same 'C astelli & Kuruczl 
(120041) atmosphere that was shown in Figure [T] As we see 
in Figure [16] the model provides a reasonable match to many 
of the effects of metallicity on the atmosphere by changing 
the distance/radius, temperature and reddening of the star. 
In particular, the metal-line blanketing effects blueward of 
^ 0.45 /im and the molecular bands in the IR are fit reason- 
ably well. The introduced shifts in the model parameters are 
(A//,A£(B-y),Af) = (-0.152,-0.166,-0.053) per dex in 
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Figure 17. Theoretical residuals from a fit to the magnitude difference due 
to metallicity between model atmospheres with A[M/H] = 0.5 dex (Am = 
'''[M/H]=-0.5 ~'"[M/H]=o.o) in sach of our filters with respect to a model that al- 
lows only for changes in distance, reddening and temperature of the Cepheid 
(Eq. 1231 ). Here we convolved Am, and /3 with a tophat function with the 
central wavelength and width of each filter, and performed the fit weighting 
each filter by the number of measurements in our dataset. 
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Figure 18. Top: Actual mean fit residuals as a function of wavelength for the 
Galaxy (red open squares), LMC (green filled circles), and SMC (blue open 
triangles). The vertical error bars are the uncertainties in the mean. Bottom: 
Relative differences of the residuals of the LMC and SMC with respect to the 
Galaxy. 



the sense that at fixed L and T a metal-poor atmosphere can 
be mimicked by making the Cepheid fainter (more distant or 
smaller), hotter, and more reddened than the same Cepheid 
with higher metallicity. These changes are then superposed 
on any shifts in the period-radius and period-temperature re- 
lations. A formalism related to Equation ( l23T l can be extended 
to determine the optimal set of filters for determinations of the 
distance, reddening, temperature and metallicity, as outlined 
in Appendix lAl 

4.2. Metallicity Effects in the Fit Residuals 



We start looking for signs of composition effects in the fit 
residuals, but it is again important to remember that most of 
the effects of metallicity are automatically projected onto the 
individual parameters of the Cepheids (f, p, E(B-V), fi). We 
discuss m etalli city-dependent signatures in these parameters 
in Section 14.41 The fit residuals contain only the manifesta- 
tions of metallicity that are orthogonal to changes in distance, 
radius, reddening and temperature. They can be identified 
by examining the residuals of the fits as a function of wave- 
length and metallicity. In order to get an idea of the expected 
signal, we repeated the procedure from Section 14.11 but in- 
stead of integrating over the model atmosphere, we convolved 
the m[ivi/Hi=-n.s -w^tM/H bo.n "sp ectra" and (3 constructed from 
the Castelli & Kurucg ( 120041) model atmospheres, and the 
ICardelh et al.. (il989>) = 3.3 extinction law with tophat 
functions having the central wavelength and width of each of 
our filters, and performed the fit in Equation (|23] | for these 
synthetic bandpasses. Individual filters were weighted by the 
number of measurements for each filter in our dataset (see 
Table We did not include the Hipparcos and Tycho fil- 
ters as there are no LMC or SMC data in these passbands. 
The residuals of this fit as a function wavelength are shown 
in Figure [TT] With these changes in the weight at any given 
wavelength, the shifts in stellar parameters are (A/i, A£'(B- 
y), At) = (-0.199, -0.208, -0.064)dex-', broadly consistent 
with the uniformly weighted theoretical comparison in Sec- 
tion l4. 1 1 We see that changing the metallicity of the stellar at- 
mosphere by 0.5 dex is largely compensated for by changes in 
the distance, reddening and temperature, leading to residuals 
of < 0.05 mag in the majority of the filters. That the residu- 
als are so small despite fitting 29 bands from 0.3 /im to 8 fim 
helps to explain why pinning down the effects of composition 
is so difficult! 

In Figure [18] we show the actual fit residuals as a function 
of wavelength A for each of the three galaxies. In order to fa- 
cilitate comparison to our theoretical model in Figure [17] the 
bottom panel shows differences in the residuals of the LMC 
and SMC from the Galaxy. We see that our results are quali- 
tatively similar to Figure [17] In particular, for »jsmc - wJCaiaxy 
we see that for A < 0.5 fim the residuals are negative and de- 
creasing with decreasing A, are approximately constant and 
positive between 0.6 /im and 1.0/im, followed by negative 
residuals in the H and K bands. The data for the LMC do 
not show such a clear trend, but the LMC metallicity is closer 
to the Galaxy and should show smaller differences than the 
SMC. Although the residuals here differ in detail from Fig- 
ure [TT] (because of the mismatch between theoretical and ac- 
tual profiles of ^ and /3 (Figures [7] and [8]l and weighting by 
the data), we see similar patterns. In particular, when compar- 
ing metal-poor to metal-rich we see negative residuals blue- 
wards of ^ 0.45 /im, slightly negative and flat residuals be- 
tween 0.45 /im and 1 .0 /im, and negative residuals in the near- 
IR 

IScowcroft et al.l(l201 Ih suggest that differences between the 
Spitzer [3.6] and [4.5] bands are a potential metallicity indi- 
cator because the [4.5] band includes the CO bandhead near 
5 /im. In Figure [18] we see that there are indeed substantial 
differences between the Galaxy and the LMC and SMC. The 
differences are smallest for the [3.6] band and increase for the 
[4.5], [5.8], [8.0] bands. The pattern of these differences does 
not agree with what we predicted in Figure[T7] — not only is the 
offset larger but the sign has rever sed. This could be a nother 
symptom of the problems in the ICastelh & Kuruc d (120041) 
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Figure 19. Light curves of two Galactic Cepheids Z Lac and VX Per with 
essentially identical periods, AP/P~ 8 X 10"^ an d metalhcities ([Fe/H] = 
0.01 and -0.05 for Z Lac and VX Per, respectively. lAndrievskv et alJ2002D . 
The light curve shapes differ due to the well-known resonance between the 
fundamental mode and the second overtone. 



model molecular opacities we discussed in Section 13.21 but 
it could also be created by different calibrations of the IRAC 
data, as the data for each galaxy are from different groups 
of authors. Nonetheless, we conclude that the pattern of the 
residuals is likely evidence for composition effects. Quantita- 
tively interpreting Figure[T8]is non-trivial, not because of any 
uncertainties in the model, but because of the effects created 
by the different weightings of the filters in the Galactic, LMC, 
and SMC data sets. 

Finally, we return to the phase dependent structures visible 
in Figure[3]that we noted in Section [33] Some of the residuals 
are due to limitations in the underlying ansatz that the light 
curve structure is solely determined by P, A^, f and p (see 
Section 12. Il l, but some of the residuals are caused by physi- 
cal differences between otherwise almost identical Cepheids. 
Figure [19] shows the light curves of two Galactic Cepheids 
whose relative period difference is only AP/P = 8 x 10"^. 
There is a significant differen ce between phases 0.5 a nd 0.7, 
where Z Lac ([Fe/H] = 0.01. lAndrievskv et alj|2002b shows 
a bump, while VX Per ([Fe/H] = -0.05) exhibits a dip. This 
is a demonstration of the well-known strong resonance be- 
tween the fundam ental mode and second overtone for 10 < 
P < 13 days (e.g. ISimon & Led 119811: lAntonello & Morellil 
[19961) . This resonance is responsible for most of the struc- 
ture seen in the residuals in Figure [3] and for the "jump" in 
the upper envelope of amplitudes at ^ 13 days in Figure [15] 
Similar (but smaller) discrepancies, especially in the depth of 
the dip preceding the rise to the maximum are also seen for 
longer periods, but we could not find any discernible pattern. 
While these resonance effects introduce noise, they should not 
significantly affect the mean properties of the Cepheids or the 
global variables as the fit residuals are quite small. We thus 
do not discuss this issue in more detail. 

4.3. Metallicity Effects on the Global Parameters 

Additional signs of the effects of composition can be 
searched for in the global parameters of our fit. Metallicity 
effects can modify the PL relations of the individual galax- 
ies, where the zero point comparison is degenerate with dis- 
tance uncertainties, but the slopes are not. The PL relations 
found by Unear least squares fits to the mean extinction- and 
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^.ML = (-4.53 ± 0. 1 1 ) - (2.76 ± 0.30) log j j(24a) 

^LMC = (_4. 16 ± 0.05) - (2.70 ± 0. 15) log [^^^ (24b) 

M!;''"''^ = (-3.93±0.07)-(2.35±0.29)log( — ) j(24c) 

\10d/ 

where My is the absolute magnitude in the V band. All stars 
within each galaxy had a uniform weight in the fit and the er- 
ror estimates from the fits are very close to those obtained by 
bootstrap resampling. The uncertainties do not include covari- 
ances in our model, unlike the PL relation in Equation (l2Ti ). 
The rms scatter of the Cepheids about these mean relations are 
o-SMC ^ 0.43, a^^'^ = 0.35, and ct'^"''"'^ = 0.43 mag. However, 
we find that the exact values of the slopes and to a lesser extent 
the zero points depend sensitively on the statistical methods 
used for performing the linear regression. 

There is a significant difference between the zero points in 
the sense that SMC Cepheids are brighter than LMC Cepheids 
which are in turn brighter than Cepheids in the Galaxy. This 
is in line with theoretical inferences from atmosphere and 
pulsational models, which predict that metal poor Cepheids 
are brighter. Arrows in Figure [14] indicate the direction of 
zero-point shifts in the V-ba nd PL relation due to metallicity 
as predicted in Section 14.11 We see that the implied shifts 
are compatible with what is observed. This issue is dis- 
cussed in greater detail in Section l4~4l and in Figure [20] An- 
other issue to consider for differences between the galaxies 
are our various distance priors. We find that the mean dis- 
tance modulus of SMC Cepheids is ^smc = 18.95 ± 0.02 mag 
with a scatter of 0.11 mag, which is very close to our prior 
of /isMC = 18.90mag with a width of O.lOmag. Given the 
LMC distance modulus of /j-lmc = 18. 50 mag, the difference 
is /^SMC - Mlmc = 0.45 ± 0.02 mag. This is in close agree- 
ment with 0.44 ± 0.05 mag determined bvlCioni et al J (120001) 
from a large sample of TRGB stars. I Storm et al.l ( 201 lb ) 
found /isMC - Mlmc = 0.47 ±0.15 mag using the infrared sur- 
face brightness method. For the Galactic Cepheids, where we 
did not use any distance priors, the distances are constrained 
by the implicit Baade-Wesselink aspects of our model, so it 
is interesting that we see only a marginal zero point shift be- 
tween the Galaxy and the LMC as we also saw in our ear- 
lier comparison to Storm etal] (1201 lajbl) . Our sample con- 
tains two Galactic Cepheids with parallax measurement and 
we found in Section [33] that our distances are in good agree- 
ment with parallax measurements. We see in Equations ( |24] | 
that the slopes agree within their uncertainties with some ev- 
idence that the slope of Galactic Cepheids is somewhat shal- 
lower. However, the worry here is that recovering the PL 
from Galactic Cepheids with widely varying distances and 
reddenings is quite difficult, and our sample has no Galactic 
Cepheids with P > 50 days. Adding both short and long pe- 
riod stars would help us to better characterize the differences 
in slopes. In summary, there are hints of differences in the PL 
relation slopes, but this is not the best probe for composition 
effects given the nature of the data. 

4.4. Metallicity Effects on the Individual Cepheid Parameters 

As was shown in Section l4~2l most metallicity effects on 
the stellar atmospheres will be absorbed into changes in dis- 
tance/radius, temperature and extinction. There are also the 
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Figure 20. Residuals with respect to the radius and temperature priors. The meaning of the symbols is the same as in Figure [T2I The individual priors and their 
widths are shown with grey dashed lines and grey dotted lines, respectively. The histograms on the bottom and left sides of the plot show distributions of p-(p) 
and f - (f ) for each of the galaxies normahzed by the number of Cepheids in each galaxy. The grey arrows in the lower left corner indicate the directions in 
which changes in p and f exactly cancel each other and there is no net change of m, for the U, V, and [3.6] bands as well as for the bolometric luminosity. The 
arrows above the temperature histogram show the expected change of f and p as a result of a metallicity increase of 0.5 dex due to changes in the atmospheric 
metallicity at fixed temperature (upper arrow) and changes in the physical properties of the star at fixed period based on pulsational models (lower arrow). 



direct changes in the mean T and R at fixed P, which also 
cause shifts in E{B-V) and jj. when a single PL relation 
({p(P)) and (f(f ))) is used. Without independent extinction 
scales, we have no means of confirming composition effects 
on extinction inferences given the very different intrinsic ex- 
tinctions of the three galaxies seen in Figure [TO] In Fig- 
ure |20] we show the residuals of the Cepheid radii and tem- 
peratures from their respective updated priors, p - (p) and 
f - (r) (Eq.l20]i. In this way, most of the period dependence 
of radius and temperature are removed and we can search 
for systematic differences between the individual galaxies. 
As a reminder, based on the discussions in Sections 14.11 and 
1421 we expect shifts of order (Ap, Af) ~ (0.04, -0.06) dex"' 
from the changes in the stellar atmospheres and ( Ap, Af) ~ 
(-0.03, -0.02) dex"' from the shifts in the pulsational models. 

Looking at the histogram of the radii, we see that SMC 
Cepheid radii have large errors, while the LMC and Galac- 
tic Cepheids generally have much smaller errors. There is no 
noticeable shift in mean radii between the individual galaxies. 
Unfortunately, the three galaxies have very different distance 
priors and available numbers of radial velocities. Since radius 
estimates will be strongly correlated with distance, the differ- 
ences in the radius residuals are dominated by these system- 
atic issues. This is particularly visible for the SMC, where 
the lack of radial velocity measurements means that the ra- 
dius determinations are dominated by the prior on p and the 
Gaussian prior on the SMC distance. This explains the rela- 



tively large errors on the SMC Cepheid radii. In the LMC, the 
distances are fixed and any variations in the mean luminosity 
are thus absorbed by the radius. Hence we get small errors 
in the radius as the mean luminosity is well determined and 
any variations in luminosity are absorbed into a tightly con- 
strained radius. In essence, these three galaxies are not ideal 
for examining metallicity effects on radius. Moreover, we see 
that the combined effects of shifts in the period-radius rela- 
tion and the changes in the radius estimate created by changes 
in the stellar atmospheres tend to cancel. 

The apparent temperature distributions are a more promis- 
ing area to look for differences because the broad wave- 
length baselines and well constrained extinctions lead to well- 
constrained temperature estimates up to the shifts created by 
changes in the stellar atmospheres (Fig.fTOli. While the distri- 
butions are broad, we see a clear trend that Galactic Cepheids 
are generally cooler than the LMC Cepheids and the SMC 
Cepheids are generally hotter. Errors on the temperatures are 
comparable for all galaxies and are smaller than the width of 
the temperature prior These relative shifts in temperature are 
robust and do not depend on the choice of the priors that fix 
the temperature scale. The simplest explanation of this pattern 
is the effect of composition as outlined in Section l4n 

Figure |20] also shows directions in the radius-temperature 
space along which there are no changes in luminosity for sev- 
eral filters as well as for the bolometric luminosity. There 
is considerably more freedom for the fit along these direc- 
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Figure 21. Absolute values of the higher-order terms compared to the lead- 
ing temperature term 2.5/3f . The top panel compares it to the typical scale of 
the second temperature derivative F^t and the temperature/metallicity cross- 
derivative r^[]y[ynj. The middle panel compares it to the typical scale of 

the first (r[M/H]) ™d second (r[[yi/H][M/H]) derivatives of M with respect to 
metalHcity [M/H]. The bottom panel compares 2.5/3-f to the derivative Tiogg 
with respect to surface gravity log g and to the change of the reddening coef- 
ficients in filter i with temperature T^-. The grey vertical bands show the 
pass bands of several common filters. 

tions because the luminosities are very well determined. We 
also show the shifts in temperature due to changes in the 
stellar temperature at fixed period and from the changes in 
the model atmosphere at fixed temperature. These two ef- 
fects have the same sign, leading to a net effect of order 
At ~ 0.12 dex"'. This broadly agrees with the observed 
shifts, given that the metallicity difference between the SMC 
and the LMC is A[Fe/H] ~ 0.4 and between the LM C and 
the Galaxy is A[Fe/H] - 0.3 dRomaniello et al.l l2008). 

Finally, it is known that the light curve amplitudes of 
SMC Cepheids are smaller than those of their LMC counter- 
parts and that this effect might be related to the metallicity 
differen ces (van Genderen 1978; Paczvnski & Pindor 2000t 
iKlagyivik & Szabadosi i2009i: iSzabados & Klagvivik i201 ll) . 
This is apparent in Figure [15] which confirms that SMC 
Cepheids with periods between 13 and 30 days do not reach 
amplitudes as high as in the LMC and Galaxy. There does 
not seem to be any discernible difference in amplitudes be- 
tween the LMC and the Galaxy. Unfortunately, the scatter 
is such that the statistical pattern is not useful for individual 
Cepheids. 

4.5. The Importance of Higher- Order Terms 

The results of Sections l4n44.4l stronglv suggest that adding 
a metallicity term to our model is a next logical step. We 
should, however, consider whether any other higher-order 
terms are equally important. These can come in the form of 
additional terms in Equation ([T]i or in more complex period- 
radius and period-temperature priors (Eq. ||6l). Here we fo- 
cus on the first question. The typical value of the first-order 



metallicity correction, essentially 9M/9[M/H], is 
d\ogF 



[M/H] ■ 



-2.5 



9[M/H] 



A[M/H] = 7A[M/H], (25) 



where 7 = -2.591ogF/9[M/H] is the "standard" metalUcity 
correction vector and A[M/H] = 0.5 dex is the typical spread 
in metallicity of Cepheid samples. Another missing first-order 
term is the surface gravity correction 
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(26) 



where Alogg = 0.5 is the typical spread in Cepheid surface 
gravities. Possible second-order terms in temperature and 
metallicity are the second derivative of the spectrum with re- 
spect to temperature (i.e. the term proportional to d(3/df) 



la^iogF 
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(27) 



the metallicity dependence of the temperature vector (the term 
proportional to d/3/d[M/n]) 



f[M/H] ■ 
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(28) 



and the second derivative of the spectrum with respect to 
metallicity 



r[M/H][M/H] - ~2.5 



1 d^logF 

2 5[M/H]^ 



(A[M/U]f 



(29) 



Here, r = 0. 1 is the typical spread in Cepheid temperatures for 
10<P< 100 days (Fig.[T3l). 

Finally, the reddening coefficients will also change with the 
temperature of the star, because the flux and the correspond- 
ing wavelength- weighted reddening will shift within the filter 
pass band. We thus consider a reddening-temperature cross- 
term 

^^'^%AE(B-V), (30) 



r^r - ■ 



df 



where AE{B-V) = 0.2 mag is the assumed spread in redden- 
ings, and ^i{f) is defined as 



-2.5 



AE(B-V) 



log 



7Q°°A-'n,(A)F.(A,f)g-^MA' 
/„°°A-in,(A)F,(A,f)dA 



(31) 

Here, n,(A) is the transmission function of filter /, which we 
assume to be a top hat for simplicity, /v is the energy flux per 
unit frequency, and is the optical depth to dust at wave- 
length A defined as £gxAE(B-V) = -2.51oge-^\ 

In Figure |2T| we compare these terms to the first-order tem- 
perature term -2.5/3f . Here, we obtained our results by eval- 
uating Equ ations (l25])-(l30b us i ng the theoretical model atmo- 
spheres of lCastelli &Kuruczl ( 12004 ) with T = 5250 KQ and 
log^ = 1.5. The first point to note is that the first order 
temperature term we use for our models is clearly the dom- 
inant term except at very short A, which is why our simple 
ansatz works so well. We see that the dominant higher-order 
term is the second derivative with respect to temperature Tff, 

^ Here we choose a different value for the temperature than in the rest 
of the paper {Tp = 5400 K) in order to avoid interpolating in the grid of 
ICastelli & Kuruc j (20041) model atmospheres. We would get essentially the 
same results with calculation centered on Tq = 5400 K. 
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which is a 10-20% correction over most of the wavelength 
range considered. Its importance grows in the near-UV, where 
it can be almost as important as the first-order temperature 
term. The metallicity dependence of the zero point FfM/H] and 
the metallicity-temperature cross-term r^[|yi/n] are modestly 
smaller, typically being ^3% corrections except in the UV. 
The second metallicity derivative r[M/H][M/H]^ the reddening- 
temperature cross-term T^f and the zero point dependence on 
surface gravity Tiogg are generally small but represent ^1% 
corrections in the optical bands. It is interesting to note that 
some of these terms may also contribute to variance in the 
light curve shapes, since any higher order terms involving the 
temperature or first order terms in log^ must also modify the 
light curves. 

5. CONCLUSIONS 

In this work, we present a method for self-consistently de- 
termining distances, reddenings, radii, and temperatures of in- 
dividual Cepheids along with the mean reddening law and 
phase-dependent variations of radius and temperature as a 
function of period. Our approach is in some senses a global 
version of the Baade-Wesselink method and it is in some 
sense s an implementation of the ideas of Freedman & Madore 
(|2010b). Our method provides statistically consistent error 
estimates for all parameters of the model (Section |2]i. We fit 
our physical model to ~5,000 radial velocity and ~ 177,000 
magnitude measurements in 29 bandpasses covering wave- 
length range from 0.3 /im to 8 /im for 287 Cepheids from the 
Galaxy, LMC, and SMC. With only four to six variables per 
Cepheid, depending on the availability of velocity data and 
the amount of freedom in the distance estimates, we fit the 
complete phase dependent data set with a magnitude scatter of 
0.051 mag and radial velocity scatter of 3.5 km s"', as shown 
in Figures |3] and m 

Our templates are built to span the period range 10 < P < 
100 days and can be used to produce model light curves or PL 
relations for any of the 29 bands. Predicting the light curve 
for an uncalibrated band is straightforward provided that the 
value of a single parameter the logarithmic flux derivative 
with respect to the temperature (Eq. |[3l. Figure |7]i, can be 
calculated or estimated for the new filter. The templates do not 
provide perfect fits, particularly for periods between 10 and 
13 days where strong resonance effects are present. The mean 
phase radius and temperature variations we obtain match the 
Hertzsprung progression in the longer period stars. 

We obtain a spectral energy distribution (Figure |6j and the 
logarithmic flux derivative with respect to the temperature /3, 
(Figure 17) for a "mean" Cepheid. When we compare our re- 
sults with estimates based on the black body law or theoretical 
stellar atmosphere models, we completely rule out black body 
models. We find better agreement with the'Castelli & Kurucz' 
(|2()04) model atmospheres, but also strong evidence for prob- 
lems in the atm osphere models (see Fremaux et al.l I2006i: 
iHauschildt et al.>,1999a.b ). The coefficients which corre- 
spond to the logarithmic derivative of the spectrum with re- 
spect to temperature, deviate from the black body law due to 
metal line blanketing and molecular opacities again in gen- 
eral agreement with theoretical atmosphere models. We are 
even able to resolve several spectral features. However, the 
theoretical atmosphere models are inadequate in describing 
the infrared spectrum, probably indicating problems with the 
molecular opacities in the models. There are no difficulties 
separating temperature and reddening in our models (Fig- 
ures |9] and [TO]i up to degeneracies created by the effect of 



composition. 

We obtain a very high precision measurement of the mean 
extinction law for 29 filters from 0.3 /um to 8 /um. The mean 
reddening law of Cepheids i n our sample (Figure [8]l shows 
departures from the .^y = 3.3 iCardelli et aU (11989 ') law usu- 
ally assumed for Cepheids. In part this is due to averag- 
ing over the three galaxies, but it clearly illustrates that ex- 
tinction laws need to be considered in efforts to increase 
the precision of Cepheid distances. This is true despite 
the trend towards infrared wavelengths, because the near- 
IR ^ values have larger fractional uncertainties. Our re- 
sults for the extinction curve also let us directly evaluate 
the Wessenheit factor TZvir = I {^%v-^ir) commonly used 
in Cepheid studies (e.g. [M adore 1982; Mad ore & FreedmanI 
[19911 120091: iNeeow & Kanbur 2005). For = 3.3 (3.1) we 
find fliat Uvic = 2.521 ±0.005 (2.475 ±0.005), where the er- 
ror comes only from the determination of ^/p, because we 
hold My fixed. The result for = 3.3 is significa ntly differ- 
ent from the usually a ssumed 2.45 iWozniak & Sta nek 199^ 
iFreedman et alj[2001l) and introduces a distance shift of about 
A/i ~ 0.01 mag for a Cepheid with E{V -Iq) = O.lOmag. 
More generally, we can calculate the Wessenheit factors for 
any combination of filters and we can correctly include the 
mutual covari ances between the fi lters in the error estimate. 
For example, iRiess et al] ( 1201 lah assumed the optical/near- 
IR Wessenheit factor of Hhvic = ^hK^v-^Ic) = 0.410 to 
reduce the uncertainty on the measurement of the Hubble 
constant. From our analysis, we find TZhvic = 0.322 ± 0.006 
(0.290 ±0.006) again at fixed ."^y = 3.3 (3 1), inconsistent 
with the assumptions by Riess et al] (1201 lal) in either case. 
This illustrates that extinction laws need to be better under- 
stood to correctly measure precise distances in the universe, 
particularly since these estimates still all assume = 3.3 (or 
3.1) and so underestimate the uncertainties. 

We find weak evidence of metallicity effects in the fit resid- 
uals. They are weak because much of the metallicity effect 
can be mimicked by shifts in the distance, reddening, and 
temperature of a Cepheid (Fig. [T6] l. The effect of metallic- 
ity on stellar atmospheres that is orthogonal to these param- 
eters only leads to residuals < 0.05 mag for a metallicity dif- 
ference of 0.5 dex (Fig.[T7l) which is approaching the regime 
where we may also need to consider absolute calibration dif- 
ferences arising from the heterogeneous data sets. The clear- 
est effect of metallicity is seen in the temperature distribu- 
tions of the Cepheids, because with 29 bands the temperature 
is well constrained independently of the extinction and dis- 
tance/radius. The SMC Cepheid temperatures are typically 
higher than those of the LMC, which are in turn higher than 
for the Galaxy, although the distributions overlap. This can 
be explained as a combination of projecting the metallicity 
effects on stellar atmospheres into changes in temperature as 
well as shifts in the period-temperature relation. The shifts 
we observe are roughly consistent with expectations. There 
may also be small shifts in radius (Fig.[20ll, but for our sample 
this is difficult to disentangle from the effects of distance pri- 
ors and would be better addressed given samples of Cepheids 
truly at a common distance and with a range of metallicities. 
Similarly, we find small differences in the zero-points of our 
PL relations, which can also be attributed to the projection of 
metallicity effects, but this has many of the same limitations 
as the search for shifts in radius. We do find shifts in the PL 
relation slopes, but our limited period range (P > 10 days) is 
not well suited for a robust investigation of slopes. 

In Appendix [A| we outHne a procedure to select optimal 
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filter sets for disentangling distance, temperature, reddening 
and metallicity or to maximize or minimize metallicity effects 
when one is concerned solely with distance determination. 
We find that sets of filters spanning the longest wavelength 
range possible, for example f/V [4.5], yield the best fit results 
for parameter estimation. Such filter combinations produce 
parameter uncertainties 0.5 to 1 order of magnitude smaller 
than the BVIc combination originally planned for the HST 
Key Project. We find that the maximum effect on distance de- 
termination due to metallicity can be expected for filter com- 
binations like UBV and BVRq- On the other hand, the slope 
of the dependence of distance on metallicity can be reduced 
to less than 10"^ magdex"' for filter combinations like VHK 
or /c[3.6][4.5]. However, some filter combinations with in- 
frared filters are quite sensitive to metallicity (e.g. UVJ). For 
these metrics, there are no particular benefits to the mid-lR 
over the near-lR, so Cepheid studies in the JWST era might 
better focus on the near-lR where Cepheids are brighter and 
the PSF will be more compact (to reduce systematic errors 
from crowding and blending). Also, mid-lR Cepheid mea- 
surements might be affected by circum stellar dust (see, e.g. 
iKervella et al.ll2006l: lOalle nne et aL 2011 ). These statements 
assume the validity of the Castelli & Kurucz (2004) model 
atmospheres, which we did find are partly inconsistent with 
the Cepheid data, so filter choices should be evaluated for a 
broader range of atmosphere models. 

Our goal in this paper was simply to carry out a complete 
analysis with the simplest possible model. That the model 
does so well helps to explain why efforts to identify higher- 
order corrections such as metallicity are s o challenging. At 
least based on the lCastelh & Kurucd (12004ft atmosphere mod- 
els, the next most important terms are the second derivative 
of the zero point with respect to the temperature followed by 
the first-order metallicity correction and the metallicity cor- 
rection to the logarithmic flux derivative with respect to the 
temperature. The second derivative of the zero point with re- 
spect to metallicity, the first order surface gravity derivative 
and the dependence of the extinction correction on tem- 
perature are relatively unimportant unless working in the UV. 
There is also a clear metallicity dependence to the mean tem- 
perature, which also indicates the need for higher-order terms 
in the period-radius and period-temperature relations. These 
estimates of higher-order eff ects are again based on th e theo- 
retical atmosphere models of ' Castelli & Kuruczl (120041) . 

There is clearly a broad range of possible future exten- 
sions. Adding additional data or galaxies is trivial, as is ex- 
tending to shorter periods. It is also possible to fit for sig- 
natures of a companion to the Cepheid which will distort the 
SED. The virtue of our approach, particularly in a cosmolog- 
ical setting, is that it forces a correct use of prior assump- 
tions in the parameter estimation and uncertainties, allows 
the data to overrule those assumptions, if necessary, and pro- 
duces error estimates incorporating the full uncertainties of 
the entire model. Improved stellar atmosphere models would 
lead to better priors, although this is less of a concern given 
the demonstrated ability of the data to constrain the SED of 
Cepheids and their dependence on temperature. Particularly 
with better atmosphere models it should be possible to re- 
lax the assumptions about extinction (.^v = 3.3 and mean 
(E(B -V)) = 0.147 mag for the LMC based on lUdalski et"an 
Il999h and put the Cepheids on their own absolute extinction 
scale. This would particularly help in finding the expected 
environmental dependencies of the extinction law and exam- 



ining their effects on distance uncertainties. 

The final imp ortant lesson from th is study, as also previ- 
ously stated by Freedman & Mador3 (l20 10b). is that there is 
no need to obtain a light curve in order to add almost all the 
information available from observing at a new wavelength. 
Given one band with enough data to determine the period and 
phase of the Cepheid at the time of any observation at an- 
other band, a single observation in this new band suffices be- 
cause the model correctly includes all the phase correlations 
between the bands. In practice, it would be wise to obtain 
at least two observations so as to have an internal check on 
the results. This removes any major barrier to systematically 
using the short wavelength filters that provide the greatest 
leverage for the control and measurement of systematic er- 
rors. Much of the drive towards longer wavelengths has been 
that with lower amplitudes the light curves can be more poorly 
sampled to still yield an accurate mean magnitude. With our 
models, any band can be poorly sampled provided its phase is 
well-determined. In s ome sense, this is a multi-band version 
of the method used bv lGerke et al.l (1201 ll) . where Cepheids in 
M81 were identified and phased with ground based images, 
and then calibrated using single epochs of HST data com- 
bined with the Stetson ( 1996) template models. Furthermore, 
when fitting an individual Cepheid, all the bands can be si- 
multaneously co-phased by using our models and measuring 
the goodness of fit to all the data as a function of period. 
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APPENDIX 



PILTER OPTIMIZATION ROR CEPHEIDS 

We can use our approach for estimating the effect of metallicity on individual filters presented in Section|4T|to design a method 
to evaluate and optimise filter choices for different metrics such as minimizing the overall error in the physical parameters or to 

i ield the s mallest or largest possib le signal of metallicity. The method is general and related to the "dark energy figure of merit" 
Albrechtet al. 2006; Wang BOOl . 

Assume that we obtain a set of mean magnitudes of a Cepheid m, where each component of m = (mi, . . . ,mN^) is a mean 
magnitude in one of the Np distinct filters. We need to determine the distance /i, reddening E(B-V), temperature t, and metallicity 
[M/H] by fitting 

m = ^ll + ME(B-V)-2.5f3T+[M/U]'y. (Al) 

Here, 1 is a vector with all components equal to unity, ^ is the reddening vector, /3 is the logarithmic change of flux with 
temperature, and 7 is the metallicity correction vector Without radial velocities, the radius is degenerate with distance and we 
merge these two quantities into fi. I n order to keep the analysis simple, we use the = 3.3 Cardelli et al.l (fl989) reddening law 
for ^ and we calculate (3 from the lCastelli & Kuruczl (2004) model atmospheres for Tq = 5400 K, logg= 1.5, and [M/H] = 0.0. 
The metallicity vector 7 is defined in Equation (|25] | from Section l43] We then convolved the models with a simple top hat model 
for the filters to compute the values of f3, and 7 for each filter 

If we optimize the fit implied by Equation (lAlb to any data set, the best fit parameters defined by minimizing the fit statistics 
are 
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where a represents the measurement errors. C is the covariance matrix of the basis functions 1, 



-2.5/3, and 7 



C = 
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(A3) 



The covariance matrix C has rank 4 and we thus require measurements in at least four filters to uniquely determine all four 
parameters. However, we have a prior knowledge on the temperature of the Cepheid (essentially the width of the instability 
strip), which we can use to reduce the number of necessary bands to three. We assume that the instability strip has a width of 
(Tt- = 0.02 mag (see Section |273T l. and we add term 1 /cr^ to the diagonal entry of C that corresponds to the temperature, which 
then reads (2.5/3)^+ 1 /a^. We assume measurement errors of cr = 0.05 mag. 

The best choice of filters then depends on the desired metric for evaluating success. One possible metric is the overall size the 
error elhpse 



£ = 1og 



1 



detC 



1/2 



(A4) 



This metric minimizes the generalized area of the error ellipse and hence tries to obtain the best joint estimate of all four param- 
eters (/.i,E(B-V),T, [M/H]). A second possible metric is 



f'=iogicr/2, 



(A5) 



which tries to obtain the smallest error in distance given that all four parameters need to be estimated. An alternate set of metrics 
is to examine the sensitivity of distance measurements to the metallicity correction 7. We rewrite Equation jAlb as 



m-[M/H]7 = /Lil-l-^E(B-V)-2.5/3r, 



(A6) 
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and interpret this equation as fitting a "true" metallicity-independent magnitude m plus a metallicity correction proportional to 
7. The solution is 



where C is the 3x3 up per-l eft submatrix of C, including the prior contribution 1 Ig\ to the diagonal temperature term. The 
second term in Equation ( IA7I ). 



is an estimate of the sensitivity of the parameter vector {fi,E{B — V),T) to metallicity effects. In particular, the change in distance 
modulus due to the change in metallicity is proportional to 5^, which has units of mag dex~'. Thus, the metric £" 



expresses the sensitivity of the distance estimate to metallicity for any given filter set. 

In order to calculate these metrics, we investigated three groups of filters. The first group has the Johnson-Cousins filters 
UBV{RI)c, the second one has Sloan filters u' g'r'i'z' , and the third includes all these filters plus the Walraven and Washington 
photometric systems. In addition, all three groups include the JHK and Spitzer IRAC bands. We first evaluated £ for all possible 
triplets of filters within each group. Table|8]shows the nine filter triplets with the smallest £ in each group normalized as £-£bvIc 
with respect to the originally planned choice of the Hubble Key project (Freedman et al. 2001), BVIq, which has £bvic = -2.379. 
The Key Project as realized used only two bands, and thus has £kp = 00 and no significant control over systematic uncertainties 
when examining individual Cepheids. We see that in all three groups, the best choices improve over the BVIq filter set by about 
0.9 dex, or about a factor of 8. The common feature of all the optimal solutions is broad wavelength coverage from the UV 
through the IRAC bands. Near-IR bands are only slightly worse than mid-IR bands. In the third group, the best filters are clearly 
dominated by the Walraven filter set. These filters were designed to be especially sensitive to temperature and metallicity of 
early-type stars and have relatively narrow widths that make them very sensitive to temperature and metallicity variations. The 
filters of the Washington system, which were designed for G and K giants, do not stand out as clearly, probably because the width 
of these filters is larger However, the improvement from these "exotic" filter sets is at most 0.14 dex, or only a factor of 1 .4. 

The fluxes of Cepheids in the near-UV bands are generally low (see Fig. |6]l and measurements in these filters are not easily 
obtained from the ground. We repeated the calculation without filters having central wavelengths shorter than 400 nm. The 
results are given in the middle section of Table [H] The best combinations simply swap U /u' for: B/ g' and we also see a decrease 
in the overall precision of about 0.4 to 0.6 dex, which corresponds to an increase in the parameter errors by a factor of 2.5 to 4.0 
compared to using the U /u' bands. Still, the total gain in precision with respect to BVIq is still 0.5 dex, which is approximately a 
factor of 3.2. 

Minimizing £ corresponds to trying to obtain the smallest joint uncertainty in all four parameters, and the filter combinations 
that perform worst are those using only near/mid-IR filters such as [4.5][5.8][8.0], 7A'[5.8], and A'[3.6][8.0]. These filter sets 
generally give £-£bvic ~ 2, which is two orders of magnitude worse than BVIc- The reason is that with wavelength coverage 
limited to the infrared, the extinction and temperature are not well constrained. We repeated the analysis allowing for four filters 
in order to see which additional filter is the most beneficial. We see that for both Johnson-Cousins and Sloan groups, adding the H 
band or a second IRAC band gives the best results. It also shows that the best complement to the Walraven and Washington filters 
is one of the Spitzer IRAC bands. Adding a fourth filter decreases £ by about 0.2 to 0.4 dex, which corresponds to a decrease in 
the parameter errors by a factor of 1 .6 to 2.5. 

If we are uninterested in any quantity other than distance, then £' may be a better metric. In Table|9]we present the filter sets 
leading to the smallest errors in distance given that the four physical parameters must be determined independently, again for the 
same three filter groups and relative to the original Key Project choice of BVIc- Again, the final Key Project filter set of only VIq 
gives £' = 00. The best results are obtained with one optical, one near-IR and one mid-IR filter, with an overall error decrease by 
about a factor of 10 relative to BVIq- Unlike the total error, the near-UV filters are not required to get the lowest error in distance 
given three filters, but they are a good addition if four filters are available. We also see that Johnson-Cousins and Sloan groups do 
not differ significantly in their errors. Adding "exotic" filter sets like the Walraven and Washington does not noticeably improve 
the distance determination. 

Finally, in Table [TOl w e present the filter sets with the smallest and largest sensitivities of the distance determination to 
metallicity £" (Eq. IIA9I ). There are filter combinations that are very insensitive to metallicity with distance shifts less than 
10""* mag dex"'. The least-metal sensitive filter combinations usually involve one optical (Rq, Ic, or V) and two infrared filters, 
which usually include at least one of the Spitzer IRAC bands. Infrared-only filter combinations like 7//[8.0] also yield small 
sensitivity to metallicity at the price of never being able to probe additional physics if it becomes necessary because they produce 
such large £ in Table |8] Interestingly, filter combinations that include the [4.5] Spitzer band, which is positio ned on t he CO 
bandhead feature and which is claimed to have potential for measuring metallicity effect in Cepheids (iScowcroft et al.ll2()l li) . can 
still have negligible metallicity sensitivity. With only three bands, it is relatively easy to construct filter combinations that are 
orthogonal to 7 after determining fi and f. 

Conversely, filter combinations that yield largest sensitivity to metallicity involve bandpasses with small Acenter- For example, 
UBV has metallicity dependence of about 1.1 mag dex"' and commonly used set BVIc has sensitivity of about 0.23 mag dex"'. 




(A7) 




(A8) 



f" = log|<5^ 



(A9) 
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Table 1 

References For Photometric Measurements. 



Reference Filters Datapoints 



_Oosterliofl" (1960) 1,2,3 351 

Walraven et al. (1964) 23,24,25,26,27 975 

Madore (1975) 1,2,3 2430 

Pel (1976) 23,24,25,26,27 5340 

Gravzeck (1978) 1,2,3 426 

Martin & Warren (1979) 1,2,3,5 1 666 

Harris (1980. 1983) 3,14,15,16,17 3679 

Martin (1980) 2 7 

_Szabados (1981) 1,2,3 1015 

van Genderen (1983) 23,24,27 1382 

Caldwell & Coulson (1984) 2,3,4,5 2489 

Mpffett& Barnes (1984) 2,3,9,10 3894 

Welch etal. (1984) 11,12,13 446 

^ulson& Caldwell (1985) 1,2,3,4,5 3904 

Coulson etal. (1985) 1,2,3.4,5 726 

Rgedman etal. (1985) 2,3,4,5 301 

Caldwell etal. (1986) 2,3,4,5 65 

Lanev&Stobie (1986, 1992) 6,7,8 3464 

Welch etal. (1987, 1993) 2,3,11,12,13 795 

Mathewson et al. (198^ 2,3 74 

^chechter etal. (1992) 11,12,13 189 

Fernie etal. (1995) 1,2,3 129 

Schmidt etal. (1995) 3,4 44 

2ebo & Wood (1995) 3,5 214 

_Bffl-nes etal. (1997) 2,3,9,10,11,12,13 641 

Perrvman&ESA(1997) 18,19,20 5954 

Kiss (1998) 1,2,3 270 

Moffettetal. (1998) 2,3,4,5 1879 

Berdnikov (1999) 1,2,3,4,5 24281 

Persson et al. (2004) 11,12, 13 4300 

Freedman et al. (2008) 33, 34, 35, 36 416 

Soszvhski et al. (2008. 2010) 3, 5 105402 

Marengo et al. (2010)" 33, 34, 35, 36 32 

N geow & Kanbur (2010) 33,34,35,36 134 

"total 177,326 



The VIqH filter combination used by lRiess et alJ (l2009aLl201 lallbh has sensitivity of about 0.006 mag dex As expected, includ- 
ing "exotic" filters like the Wakaven system greatly increases metalUcity sensitivity as these filters were designed to do exactly 
that. The original Hubble Key Project choice of BVIc has a moderate sensitivity of 0.23 mag dex"', while using only Vic has 
0.066 mag dex"'. 

This exploration does depend on the c hoices for /3, and p'. Here, we used the ICastelli & Kuruczl (120041) theoretical model 
atmospheres to produce /3 and 7, and the lCardelli et alJ (119891) = 3.3 model of We know these models have problems as 
we discussed in the main part of our study. If we use our empirical estimates of /3 and we do not see significant changes in the 
general trends outlined in this Appendix. However, since the primary question is the effects of 7, a well-designed survey should 
consider a range of models for 7 rather than the single case we considered here. 
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Table 2 

References For Radial Velocity Measurements. 



Reference Datapoints 



Joy (1937) 25 

Stibbs (1955) 62 

Feast (1967) 58 

Lloyd Evans (1968. 1980) 72 

Grayzeck(1978) 53 

Coulson & Caldwell (1985) 794 

Coulson et al. ( 1 98T) 232 

Tmbertetal. (1985) 258 

Beavers & Bitter ( 1 986) 21 

Caldwell etal. (1986) 62 

Barnes et al. (1987,_1988) 370 

Mathewson et al. (1988) 36 

Tmbertetal. (1989) 81 

Wilson etal. (1989) 63 

Metzger et al. (1991. 1992) 145 

Gorynya et al. (1992. 1996. 1998) 882 

Bersier etal. (1994) 681 

Pont etal. (1994, 1997) 189 

Szabados & Pont ( 1 998) 33 

Berdnikov (1999) 382 

Tmbert (1999) 266 

Kiss & Vinko (2000) 33 

Storm et al. (2004) 163 

^roenewegen (2008) 13 

Nardetto et al. (2009) 57 

total 57nT 



Table 3 

Filters and Filter Properties. 



No. 


Name 


Acenter [nm] 


A A [nm] 


Datapoints 


M, 


A- 




1 


V 


365 


35 


2891 


1.6600 ±0.0368 


10. 188 ±0.006 


5.246 ±0.003 


2 


B 


441 


49 


11574 


1.0306 ±0.0264 


7.329 ±0.003 


= 4.303 


3 


V 


549 


43 


25893 


0.1557±0.0186 


= 5.1404 


= 3.303 


4 


Rc 


659 


75 


7833 


-0.2792 ±0.0154 


4.229 ±0.002 


2.623 ±0.002 


5 


Ic 


806 


54 


100369 


-0.6877 ±0.0125 


3.393 ±0.002 


1.993 ±0.003 


6 


■'SAAO 


1216 


202 


1156 


-1.1488 ±0.0095 


2.422 ± 0.006 


0.784 ± 0.007 


7 


^'SAAO 


1629 


229 


1156 


-1.5679 ±0.0075 


1.719±0.0()6 


0.422 ± 0.007 


8 


^SAAO 


2205 


269 


1152 


-1.6356 ±0.0071 


1.548 ±0.006 


0.200 ±0.008 


9 


Rl 


693 


149 


1042 


-0.4787 ±0.0138 


3.750 ±0.007 


2.308 ±0.003 


10 


h 


879 


156 


1040 


-0.9137±0.0110 


2.904 ±0.007 


1.478 ±0.005 


11 


■Icno 


1246 


120 


2013 


-1.1 898 ±0.0089 


2.246 ± 0.005 


0.724 ±0.007 


12 


Hcno 


1623 


150 


2013 


-1.5571 ±0.0074 


1.724 ±0.005 


0.426 ± 0.007 


13 


Kcno 


2217 


200 


2005 


-1.6435 ±0.0070 


1.553 ±0.005 


0.255 ±0.008 


14 


C 


391 


78 


736 


1.3852±0.0316 


8.746 ±0.008 


4.725 ±0.003 


15 


M 


509 


74 


736 


0.3660 ± 0.0205 


5.647 ±0.007 


3.633 ±0.002 


16 


Ti 


633 


56 


736 


-0.2774 ±0.0152 


4. 155 ±0.007 


2.682 ±0.002 


17 


Ti 


789 


99 


736 


-0.6781 ±0.0124 


3.356 ±0.006 


1.934 ±0.003 


18 


Hp 


504 


157 


5816 


0.3534 ±0.0184 


5.080 ±0.003 


3.154±0.001 


19 


Bj 


419 


51 


69 


1.7390 ± 0.0249 


6.634 ±0.034 


3.523 ±0.018 


20 


Vt 


523 


67 


69 


0.2065 ±0.0196 


5.018 ±0.034 


3.520±0.018 


23 


Vw 


547 


60 


2268 


-6.6645 ±0.0190 


5.243 ±0.004 


3.417 ±0.001 


24 


Bw 


433 


29 


2199 


-5.7640 ±0.0286 


7.920 ±0.005 


4.5 15 ±0.002 


25 


t/w 


363 


18 


1244 


-4.3928 ±0.0335 


9.261 ±0.007 


5.254 ±0.003 


26 


Ww 


326 


10 


524 


-3.5206 ±0.0350 


9.646 ±0.011 


5.753 ±0.008 


27 


Lw 


384 


10 


1462 


-4.7697 ±0.0377 


10.402 ±0.008 


4.887 ±0.006 


33 


IRAC 3.6 ^tm 


3550 


370 


142 


-1.7672 ±0.0075 


1.442 ±0.023 


0.171 ±0.008 


34 


IRAC 4.5 ^tm 


4493 


510 


149 


-1.7263±0.0081 


1.726 ±0.020 


0.129 ±0.006 


35 


IRAC 5.8 ^m 


5731 


705 


142 


-1.7603 ±0.0082 


1.657 ±0.027 


0.086 ±0.004 


36 


IRAC 8.0 /.tm 


7782 


1440 


149 


-1.7901 ±0.0071 


1.397±0.018 


0.052 ±0.002 



Note. — List of filters along with their central wavelength Acciner, zero point M,, temperature coefficient and extinction coefficient 
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Table 4 

Template Coefficients. 



log (P/lOd) 


;■ 


Ct.; 














0.00 


1 


-0.655942 




-0.799722 




1.00000 




0.000000 




2 


0.432002E- 


■01 


-0.940700E- 


■02 


0.424302E- 


■01 


0.816138E-01 




3 


0.166285E- 


■01 


0.658372E- 


■01 


-0.105897E^ 


■01 


-0.390460E-01 




4 


0.395908E- 


■02 


-0.2888 1 3E- 


■01 


0.822020E- 


■02 


0.118774E-01 




5 


-0.588958E- 


■02 


0.805280E- 


■02 


-0.352807E^ 


■03 


-0.958669E-03 




6 


0.73575 IE- 


■02 


0.151437E- 


■02 


0.511735E- 


■02 


-0.142921E-02 




7 


-0.615025E- 


■02 


-0.39449 IE- 


■02 


-0.172599E^ 


■05 


0.27047 lE-02 




8 


0.418875E- 


■02 


0.347581E- 


■02 


0.957207E- 


■03 


-0.494847E-02 




9 


0.190958E- 


■03 


0.254600E- 


■04 


0.328564E- 


■02 


-0.339728E-04 




10 


0.704976E- 


■03 


0.614612E- 


■03 


0.373839E- 


■03 


-0.345807E-02 




11 


-0.463388E- 


■03 


-0.341031E- 


■03 


0.287978E- 


■02 


-0.331 189E-03 




12 


-0.2861 17E- 


■03 


0.178882E- 


■03 


-0.168316E^ 


■03 


0.228903E-03 




13 


-0.190860E- 


■03 


-0.288726E- 


■03 


-0.493956E^ 


■03 


0.255491E-03 




14 


-0.113362E- 


■03 


-0.314758E- 


■04 


-0.389795E^ 


■03 


-0.563042E-03 




15 


-0.254223E- 


■03 


-0.177749E- 


■03 


0.117149E- 


■02 


0.560073E-03 




16 


-0.560226E- 


■04 


0.128664E- 


■04 


-0.120511E^ 


■02 


0.640936E-03 




17 


-0.186283E- 


■03 


-0.141655E- 


■03 


-0.483337E^ 


■03 


-0.642305E-03 




18 


-0.158575E- 


■03 


-0.150494E- 


■04 


0.965609E- 


■03 


0.200574E-03 




19 


-0.768975E- 


■04 


-0.552935E- 


■04 


-0.662260E^ 


■04 


0.343053E-03 




20 


-0.230873E- 


■03 


-0.156252E- 


■03 


0.51041 IE- 


■03 


-0.5531 18E-03 


0.02 


1 


-0.654341 




-0.763829 




1.00000 




0.000000 




2 


0.238339E- 


■01 


-0.128259E- 


■01 


0.964536E- 


■02 


0.115707 




3 


0.490282E- 


■01 


0.632282E- 


■01 


-0.141993E^ 


■01 


-0.394291E-01 




4 


-0.146646E- 


■01 


-0.365181E- 


■01 


0.126668E- 


■01 


0.169339E-01 




5 


-0.527408E- 


■02 


0.141425E- 


■01 


-0.464929E^ 


■02 


-0.613523E-03 



Note. — The table lists Fourier coefficients of the templates at Np = 19 anchor points from log{P/10d) = 
0.00 to 1 .00. Each anchor point contains four columns, c-rj, s-rj, Cpj, and Spj, where each line gives the j-th 
coefficient, j = 1 . . .N^, where = 20 (see Eqs. (Tj-d))- The template at any period is obtained by linear 
interpolation of all coefficients in logP bracketing the desired period. The complete Table is provided in the 
electronic edition. 



Table 5 

Galactic Cepheids 



Name 


P[d] 


A* [mag] 


BZ Cyg 


10.1420 


11.376 ±0.040 


SY Aur 


10.1446 


11.708 ±0.040 


f Gem 


10.1505 


7.892 ±0.025 


FQCar 


10.2742 


13.236 ±0.079 


AN Aur 


10.2893 


12.981 ±0.055 


YSct 


10.3415 


11.021 ±0.039 


MZCen 


10.3538 


13.842 ±0.035 


FOCar 


10.3560 


13.582 ±0.061 


FRCar 


10.7170 


12.799 ±0.016 


TWNor 


10.7861 


11.700 ±0.032 


Z Lac 


10.8857 


11.291 ±0.100 


VX Per 


10.8866 


12.003 ±0.019 


XXCen 


10.9536 


10.669 ±0.011 


TY Set 


11.0540 


11.690 ±0.033 


SVPer 


11.1292 


12.520 ±0.023 


DR Vel 


11.2013 


11.314±0.017 


V438 Cyg 


11.2108 


11.297 ±0.039 


AD Cam 


11.2622 


13.878 ±0.042 


HZ Per 


11.2796 


13.844 ±0.033 


V340 Nor 


11.2894 


11.689 ±0.061 



£'(B-V)[mag] v[kms-l] 



X 10- 



0.916±0.008 
0.446 ±0.006 
-0.031 ±0.006 
0.597 ±0.032 
0.568 ±0.008 
0.735 ±0.005 
0.665 ± 0.007 
0.380 ±0.005 
0.285 ±0.005 
1.090 ±0.006 
0.374 ±0.003 
0.484 ±0.005 
0.238 ±0.003 
0.927 ± 0.006 
0.533 ±0.006 
0.609 ± 0.005 
1.111±0.013 
0.983 ±0.012 
1.312±0.012 
0.458 ±0.016 



-11.4±0.1 
0.0±0.1 
6.1 ±0.0 
-3.7 ±0.5 
1.9±0.2 
14.2 ±0.2 
-24.4 ±0.2 
-12.1 ±0.4 
-6.9 ±0.1 
-56.1 ±0.2 

-36.3 ±0.1 
-18.1±0.0 
25.0 ±0.2 
-5.8 ±0.1 
20.5 ±0.1 
-9.1 ±0.2 
-60.8 ±0.2 
-28.5 ±0.2 
-40.7 ±0.1 



0.782 ±0.008 
0.723 ±0.008 
0.819 ±0.005 
0.749 ±0.014 
0.814±0.011 
0.790 ± 0.008 
0.952 ±0.007 
0.846±0.012 
0.845 ±0.003 
0.797 ±0.006 
0.803 ±0.020 
0.819 ±0.004 
0.720 ±0.002 
0.791 ±0.006 
0.913 ±0.005 
0.815 ±0.003 
0.828 ±0.008 
0.826 ±0.008 
0.896 ±0.006 
0.877 ±0.012 



0.032 ±0.002 
0.049 ±0.002 
-0.003 ± 0.002 
-0.016 ±0.006 
0.036 ±0.002 
0.000 ±0.002 
-0.011 ±0.002 
0.004 ±0.002 
-0.001 ±0.002 
-0.016 ±0.002 
0.017 ±0.002 
0.026 ±0.002 
0.013 ±0.002 
0.007 ± 0.002 
0.067 ± 0.002 
-0.012 ±0.002 
-0.000 ± 0.003 
0.042 ± 0.003 
0.003 ±0.003 
0.043 ± 0.006 



2.2943 ±0.0038 
2.6711 ±0.0047 
2.2276 ±0.0036 
3.3464 ±0.0044 
2.8273 ±0.0032 
3.3299 ±0.0034 
2.8914 ±0.0047 
2.4539 ±0.0026 
3.2252 ±0.0031 
3.6645 ±0.0034 
3.7122 ±0.0028 
2.8020 ±0.0023 
3.6481 ±0.0027 
3.6861 ±0.0032 
3.0448 ±0.0029 
3.2269 ±0.0032 
3.8962 ±0.0048 
3.8923 ±0.0047 
3.4095 ±0.0044 
1.3556 ±0.0037 



Note. — For each Cepheid, the table lists the period P, distance modulus reddening EiB -V). mean radial velocity \', mean logarithmic radius p, 
mean logarithmic temperature r and the amplitude A^. The complete Table is provided in the electronic edition. 
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Table 6 

LMC Cepheids 



Name 



P [d] /t [mag] E{B - V) [mag] v [km s"' ] 



X 10^ 



HV 01016 
OIII 2562 
HV 00923 
OIII 0057 
HV 02371 
OIII 01 14 
HV 02277 
OIII 3284 
HV 05551 
HV 13064 
OIII 0016 
OIII 0978 
OIII 1784 
HV 12078 
OIII 0798 
HV 12537 
HV 02280 
HV 02432 
HV 02864 
HV 02598 
HV 00921 



10.0438 
10.0716 
10.2547 
10.3075 
10.3479 
10.3727 
10.4135 
10.4772 
10.4852 
10.4872 
10.5070 
10.5199 
10.6628 
10.6876 
10.7104 
10.7923 
10.8548 
10.9192 
10.9853 
10.9913 
11.0905 



18.480 
18.479 
18.511 
18.558 
18.514 
18.552 
18.533 
18.447 
18.544 
18.482 
18.554 
18.517 
18.492 
18.427 
18.526 
18.531 
18.525 
18.494 
18.426 
18.486 
18.506 



0.147±0.019 
0.147±0.019 
0.144 ±0.019 
0.150±0.019 
0.142±0.019 
0.142±0.019 
0.142±0.019 
0.147±0.019 
0.093 ±0.011 
0.191 ±0.019 
0.169±0.019 
0.157±0.019 
0.147±0.019 
0.088±O.01O 
0.145±0.019 
0.142±0.019 
0.143±0.019 
0.247 ±0.010 
0.095 ± 0.006 
0.138±0.019 
0.147±0.019 



259 



±0.1 



0.761 
0.788 
0.784 
0.809 
0.780 
0.789 
0.798 
0.774 
0.777 
0.841 
0.905 
0.840 
0.777 
0.802 
0.789 
0.819 
0.801 
0.787 
0.791 
0.806 
0.811 



±0.003 
±0.003 
±0.002 
±0.003 
±0.003 
± 0.003 
±0.003 
± 0.003 
±0.001 
±0.003 
±0.003 
± 0.003 
±0.003 
±0.001 
± 0.003 
±0.003 
± 0.003 
±0.002 
±0.001 
± 0.003 
±0.003 



0.025 ±0.006 
0.033 ± 0.006 
0.039 ±0.006 
0.027 ± 0.006 
0.046 ± 0.006 
0.048 ± 0.006 
0.051 ±0.006 
0.025 ±0.006 
0.016 ±0.003 
-0.099 ±0.006 
-0.004 ± 0.006 
0.013 ±0.006 
0.023 ± 0.006 
0.004 ± 0.003 
0.032 ± 0.006 
0.053 ±0.006 
0.041 ±0.006 
0.100 ±0.004 
0.020 ± 0.002 
0.060 ±0.006 
0.030 ± 0.006 



2.6563 ±0 
0.9424 ± 
2.3734 ±0 
2.0783 ±0 
3.3251 ±0 
2.5954 ±0 
3.4945 ±0 
2.3392 ±0 
2.0011 ±0 
3.0641 ±0 
0.6998 ± 0, 
1.7861 ±0, 
2.9994 ±0 
2.0235 ±0, 
1.7056 ±0, 
3.2753 ±0, 
2.9274 ±0, 
2.93 10 ±0 
3.7166 ±0 
2.1020±0, 
2. 1498 ±0, 



.0067 
.0018 
0030 
.0033 
.0042 
0040 
.0050 
0043 
.0033 
.0056 
0018 
0027 
.0042 
0030 
0023 
0037 
0032 
.0037 
.0037 
0020 
0024 



Note. 



iame as in 



TablelSl but for the LMC Cepheids. The complete Table is provided in the electronic edition. 



Table 7 

SMC Cepheids 



Name 


P[d] 




H [mag] 


E{B-V) [mag] 


V [km s ' ] 


P 


T 


A2 X 10^ 


HV 06320 


10.0998 


19.017 ±0.071 


0.148±0.015 






0.757 ±0.014 


0.078 ± 0.004 


3.7284 ±0.0095 


HV 02063 


11.1663 


IS 


.881 ±0.071 


0.1 67 ±0.007 






0.815±0.014 


0.055 ±0.002 


2.8573 ±0.0033 


HV 02017 


11.4091 


IS 


.870 ±0.071 


0.175 ±0.007 






0.823 ±0.014 


0.060 ± 0.002 


3.0375 ±0.0030 


HV 01610 


11.6431 


IS 


.944 ±0.071 


0.171 ±0.009 






0.814±0.014 


0.070 ±0.003 


3.7982 ±0.0038 


HV 00856 


12.1557 


19.038 ±0.071 


0.273 ±0.008 






0.809 ±0.014 


0.081 ±0.003 


3.6536 ±0.0043 


HV 01365 


12.4127 


19.483±0.016 


0.231 ±0.005 


126 


8±0.1 


0.885 ±0.003 


0.070 ± 0.002 


3.3213 ±0.0026 


HV 02227 


12.4663 


IS 


.866 ±0.071 


0.201 ±0.010 






0.851 ±0.014 


0.048 ± 0.003 


3.3765 ±0.0036 


HV 02230 


12.5294 


IS 


.823 ±0.071 


0.218±0.010 






0.861 ±0.014 


0.051 ±0.003 


3.0924 ±0.0067 


HV 02052 


12.5778 


IS 


.858 ±0.071 


0.084 ± 0.008 






0.855 ±0.014 


0.058 ±0.003 


3.8876 ±0.0045 


HV 01744 


12.6237 


IS 


.925 ±0.071 


0.1 35 ±0.009 


141 


9±5.8 


0.842 ±0.014 


0.061 ±0.003 


4.0334 ± 0.0046 


HV 01873 


12.9396 


19.026 ±0.071 


0.1 85 ±0.007 






0.830 ±0.014 


0.060 ± 0.002 


4.4509 ±0.0041 


HV 02225 


13.1477 


IS 


.877 ±0.071 


0.119±0.014 






0.865 ±0.014 


0.024 ± 0.004 


3.5367 ±0.0095 


HV 02202 


13.1921 


IS 


.852 ±0.071 


0.227 ±0.008 






0.871 ±0.014 


0.073 ±0.002 


3.7139 ±0.0051 


HV 00827 


13.4642 


IS 


.935 ±0.071 


0.1 86 ±0.006 






0.861 ±0.014 


0.072 ± 0.002 


3.6416 ±0.0031 


HV 02189 


13.4703 


IS 


.989 ±0.071 


0.360 ±0.022 






0.850 ±0.014 


0.1 14 ±0.007 


3.5912 ±0.0039 


HV 01326 


13.7250 


19.041 ±0.071 


0.169 ±0.009 






0.845 ±0.014 


0.049 ± 0.003 


3.4243 ±0.0033 


HV 01335 


14.3814 


19.079 ±0.071 


0.243 ±0.009 






0.852 ±0.014 


0.076 ±0.003 


3.5012 ±0.0039 


HV 02088 


14.5796 


IS 


.955 ±0.071 


0.208 ±0.008 






0.881 ±0.014 


0.052 ± 0.002 


4.4474 ± 0.0053 


HV 00843 


14.7136 


IS 


.927 ±0.071 


0.120 ±0.009 


166 


8±5.8 


0.888 ±0.014 


0.002 ± 0.003 


4. 1825 ±0.0037 


HV 02233 


15.1677 


IS 


.807 ±0.071 


0.271 ±0.016 






0.922 ±0.014 


0.102 ±0.004 


5.0158±0.0141 


HV 01442 


15.2886 


IS 


.964 ±0.071 


0.232 ±0.006 


130 


0±7.1 


0.893 ±0.014 


0.046 ± 0.002 


4.4427 ±0.0041 



Note. — Same as in Table[5] but for the SMC Cepheids. The complete Table is provided in the electronic edition. 
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Tables 

Filter combinations that yield the smallest overall error S. 



UBV(RI)cJHK+mAC 


£-£bvic 


u'g'/i'z'JHK+mAC 




All 


£--£-Bvic 


{/y[4.5] 


-0.892 


H'r'[4.5] 


-0.838 


MLw[5.8] 


-0.979 


UV[5.S] 


-0.891 


H'r'[5.8] 


-0.834 


MLw[8.0] 


-0.979 


UV[H.Q] 


-0.884 


H'r'[8.0] 


-0.823 


MLw[4.5] 


-0.977 


UV[3.6] 


-0.864 


«'g'[4.5] 


-0.811 


g'Lw[8.0] 


-0.966 


URc[4.5] 


-0.861 


M'g'[5.8] 


-0.811 


g'Lw[5.8] 


-0.966 


URc[5.S] 


-0.858 


«'g'[8.0] 


-0.806 


g'Lw[4.5] 


-0.963 


URcim 


-0.847 


M'r'[3.6] 


-0.801 


MLw[3.6] 


-0.963 


UVK 


-0.842 


u'g'[3.6] 


-0.789 


(yL)w[5.8] 


-0.959 


URcD.6] 


-0.823 


u'r'K 


-0.778 


(yL)w[4.5] 


-0.958 



BRc[4.5] 


-0.487 


g'r'[4.5] 


-0.221 


Bri[4.5] 


-0.512 


BRc[5.8] 


-0.478 


g'r'[5.8] 


-0.208 


B/[4.5] 


-0.509 


BRcim 


-0.455 


g'r'vm 


-0.179 


Bri[5.8] 


-0.505 


BV[A.5] 


-0.452 


g'r'ViM 


-0.155 


Br' [5. 8] 


-0.502 


W[5.8] 


-0.447 


/i"[4.5] 


-0.149 


Bi?c[4.5] 


-0.487 


J/y//[4.5] 


-1.057 


«'gV'[4.5] 


-1.009 


M(WL)w[8.0] 


-1.357 


J/y//[5.8] 


-1.052 


H'g'r'[5.8] 


-1.006 


M(H'L)w[5.8] 


-1.357 


J/y[4.5][8.0] 


-1.043 


«'r'ff[4.5] 


-0.998 


M(WL)w[4.5] 


-1.354 


J/y[4.5][5.8] 


-1.042 


M'g'//[4.5] 


-0.996 


g'(ll'Z.)wf8.0] 


-1.347 


f/y[3.6][4.5] 


-1.041 


M'gV'[8.0] 


-0.995 


g'(U/L)w[5.8] 


-1.346 



Note. — £ -£bv1q represents the logarithmic reduction in the parameter uncertainties compared to using the BVIc 
filters. The upper part of the Table shows results for all filters, the middle part shows results for filters with Acenter > 
0.4 jum, and the lower part for four filters in the set. 



Table 9 

Filter combinations that yield the smallest error in distance £' . 



UBV(RI)cJHK+mAC £'-£bvi u'g'/i'z'JHK+JRAC £'-£bvi All £'-£b 



VHK 


-1.145 


r'ff[3.6] 
g'ff[4.5] 
g'H[5.8] 
r'/f[8.0] 


-1.129 


VHK 


-1.145 


Rcff[3.6] 


-1.122 


-1.116 


7iff[3.6] 


-1.130 


Rcff[8.0] 


-1.104 


-1.115 


VwHK 


-1.130 


Bff[4.5] 


-1.088 


-1.109 


r'H[3.6\ 


-1.129 


/cH[8.0] 


-1.087 


/'//[3.6] 


-1.107 


RcH[3.6] 


-1.122 


/cH[3.6] 


-1.084 


i'H[S.Q] 


-1.096 


g'HH-5] 


-1.116 


BK[4.5] 


-1.031 


z'm&.O] 


-1.078 


g'H[5.8] 
r'H[&.0] 


-1.115 


RcK[&.0] 


-1.013 


z'H[3.6] 
g'KlS.S] 


-1.047 


-1.109 


RcK[3.6] 


-1.011 


-1.036 


BwH[4.5] 


-1.108 


RcHK[3.6] 


-1.211 


r'HK[3.6] 


-1.217 


(BW)v/HK 


-1.248 


UBHK 


-1.204 


g'Hm.5] 


-1.214 


(BW)w//[3.6] 


-1.245 


IcHK[3.6\ 


-1.199 


i'HK[3.6] 


-1.207 


(UB)vjHK 


-1.233 


UBH[3.6] 


-1.191 


z'HK[3.6] 
u'g'HK 


-1.197 


BW^HK 


-1.233 


VH[3.6] [5.8] 


-1.179 


-1.194 


-1.228 



Note. — -^Bvir represents the logarithmic reduction in the distance uncertainty compared to using the BVIc filters. 



The upper part of the Table shows results for all filters and the lower part for four filters in the set. 
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Table 10 

Filter combinations that yield the smallest and largest change in distance due to metallicity. 



UBV(RI)cJHK+lRAC £" u'g'r'i'z'JHK+lRAC £" All E' 



Smallest change in distance 



/c[3.6][4.5] 


-4.344 


i'H[8.0] 


-4.161 




-5.149 


VHK 


-3.677 


,-'[3.6][4.5] 


-3.970 


/c[3.6] [4.5] 
i'H[8.0] 


-4.344 


RcH[SX)] 


-3.662 


r'U [8.0] 


-3.738 


-4.161 


i?c[3.61[4.5] 


-3.599 


r'[3.6][4.5] 
z'[3.6][4.5] 


-3.409 


CWw[8.0] 


-4.057 


VJ[4.5] 


-3.371 


-3.342 


(■'[3.6] [4.5] 


-3.970 


/?c[5.8][8.0] 


-3.269 


i'[5.8][8.0] 


-3.319 


Rci'Ti 


-3.926 


/c[5.8][8.0] 


-3.255 


z'[5.8][8.0] 
r' [5.8] [8.0] 


-3.315 


T2[3.6] [4.5] 
Vr'H 


-3.840 


IcJH 


-3.215 


-3.250 


-3.836 


/ff[8.0] 


-3.148 


r'/ff 


-3.209 


r'ff[8.0] 


-3.738 


/[5.8][8.0] 


-3.129 


7H[8.0] 


-3.148 


/cM[4.5] 


-3.731 


Largest change in distance 


UBV 


0.024 


u'g'r' 


-0.094 


Uu'Uw 


0.936 


UVRc 


-0.103 


u'r'i' 


-0.301 


B{BL)w 


0.678 


UBRc 


-0.186 


u'g'i' 


-0.357 


C(WL)w 


0.632 


BVRc 


-0.327 


u'r'z' 


-0.441 


B{WL)^ 


0.508 


URqIc 


-0.338 


u'i'z' 


-0.471 


(BWL)^ 


0.453 


UVIc 


-0.384 


u'g'z' 


-0.486 


B(BW)w 
g'MLw 


0.409 


UBIc 


-0.447 


u'r'J 


-0.570 


0.363 


URcJ 


-0.564 


u'i'J 


-0.593 


UBWw 


0.360 


UVJ 


-0.575 


u'g'J 


-0.623 


UBC 


0.323 


UIcJ 


-0.581 


u'z'J 


-0.686 


BCWw 


0.314 



Note. — S" represents the logarithm of the absolute change in distance modulus with metallicity. The 
upper part of the Table shows filter sets with the smallest changes, while the lower shows those with the largest 
changes. 



